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CHAPTER  1 
GENERAL  INFORMATION 


THE  FINITE  ELEMENT  CODE  AT  I  LA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  "INST  I TUT  SUPER  I  EUR  D' ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D* ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN.  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE). 


Chapter  1.  General  information 


1.1  INTRODUCTION 


The  finite  element  code  AT  I  LA  has  been  specifically  developed 
to  aid  the  design  of  piezoelectric  devices,  mainly  for  sonar 
applications.  Thus,  It  is  able  to  perform  the  modal  analysis  of 
both  ax Isymmetr leal  and  fully  three-dimensional  piezoelectric 
transducers.  It  can  also  provide  their  harmonic  response  under 
radiating  conditions:  nearfield  and  farfield  pressure, 
transmitting  voltage  response,  directivity  pattern,  electrical 
impedance,  as  well  as  displacement  field,  nodal  plane  positions, 
stress  field  and  various  stress  criteria...  Its  accuracy  and  its 
ability  to  describe  the  physical  behavior  of  various  transducers 
(Tonpilz  transducers,  double  headmass  symmetrical  length 
expanders,  free  flooded  rings,  f lextensional  transducers,  bender 
bars,  cylindrical  and  trilaminar  hydrophones...)  have  been  checked 
by  modelling  more  than  twenty  different  structures  and  comparing 
numerical  and  experimental  results.  These  tests  were-  all 
successful,  provided  that  sufficient  care  was  previously  given  to 
the  physical  meanings  of  the  models.  Most  of  them,  as  well  as  the 
general  modelling  rules,  have  been  described  in  reports,  papers 
and  communications,  which  are  partially  listed  in  the  next 
sect  ion. 

AT  I  LA  has  been  designed  for  the  "  Groupe  d* Etudes  et  de 
Recherche  en  Detection  Sous-Marine,  G.E.R.D.S.M."  from  the 
"Direction  des  Constructions  et  Armes  Navales,  D.C.A.N.",  in 
TOULON  (FRANCE),  by  the  Acoustics  Laboratory  of  the  "Inst  I  tut 
Super leur  d'Electronique  du  Nord,  I.S.E.N.",  in  LILLE  (FRANCE).  It 
has  first  been  a  research  tool,  for  specific  optimisation 
purposes,  but,  due  to  the  request  for  Its  use  by  a  large  number  of 
people,  it  has  then  been  made  more  user  friendly  and  documented, 
while  a  specific  effort  has  been  carried  out  to  Improve  Its 
versatility  as  well  as  to  provide  more  efficient  graphic  display 
abilities.  At  the  present  time,  AT  I  LA  has  been  run  on  various  IBM, 
VAX  and  UN  I  VAC  mainframes,  as  we  I  I  as  a  FPS  array  processor  and 
VAX  and  MASSCOMP  workstations.  Nevertheless,  AT  I  LA  is  not  a  "black 
boje"  and  users  are  assumed  to  have  a  firm  background  in  finite 
element  analysis. 

ATILA  has  been  developed  over  the  past  8  years  by  a  lot  of 
people,  but  mainly  by  K.  ANIFRANI,  R.  BOSSUT,  J.L.  CARTON,  J.C. 
DEBUS,  J.N.  DECARPIGNY,  B.  DUBUS,  R.  EDDE,  B.  HAMONIC,  D.  MOREL, 
P.  TIERCE  at  I.S.E.N.,  and  B.  TOCQUET  and  D.  BOUCHER  at 
G.E.R.D.S.M..  Current  developments,  Including  new  algorithms,  new 
material  description,  material  loss  description  and  finite 
element/ integral  equations  coupling  are  proceeding  at  I.S.E.N., 
while  software  engineering  is  now  provided  by  the  SINAPTEC 
company,  LILLE  (FRANCE). 


Information  and  help 
specific  conditions,  from: 
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about  AT  I  LA  can  be  obtained,  under 


J.N.  DECARPIGNY 

Acoustics  Laboratory 

Instltut  Super leur  d'E lectronlque  du  Nord 

41  boulevard  Vauban 

59046  LILLE  Cedex .  FRANCE 

Tel:  (33). 20. 30. 62. 20    Telex:  GEFIRN  120.369 


P.  TIERCE 

SINAPTEC  Company 
41  boulevard  Vauban 
59800  LILLE.  FRANCE 
Tel:  (33). 20. 57. 80. 04 


Telex:  GEFIRN  120.369 


Following  this  section  and  the  reference  list,  the  reader  will 
find  the  description  of  the  successive  stages  of  an  AT  I  LA  job  and 
their  connection  to  the  manual  organisation.  Then,  detailed 
information  is  collected  In  five  chapters. 
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1.3  GENERAL  ORGANISATION  OF  AN  AT  I  LA  JOB 


Running  an  ATILA  Job  requires  different  well  defined  stages. 
This  section  describes  briefly  their  organisation  and  their 
connections  to  the  various  following  chapters  of  this  manual. 

1.3.1.  THE  MODEL  DEFINITION. 

Being  given  the  physical  problem,  the  user  has  to  define  the 
type  of  analysis  which  is  required:  computation  of  the  static 
deformation  of  an  elastic  stucture  under  lumped  mechanical  load, 
modal  analysis  of  a  piezoelectric  stack  with  short  circuit 
electrical  condition,  computation  of  the  transmitting  voltage 
response  of  a  transducer...  This  type  of  analysis  has  to  be 
related  to  solving  algorithms  which  are  listed  and  described  in 
chapter  2.  New  solvers  are  currently  under  test  or  development  in 
I. S.E.N.  Acoustics  laboratory  and  detailed  Information  can  be 
provided  upon  request.  The  user  has  also  to  determine  the  type(s) 
of  elements  which  is  (are)  required:  three-dimensional  elements, 
ax isymmetr leal ,  plane  stress  or  plane  strain  elements,  plate  or 
shell  elements  relying  on  specific  approximations...  The  ATILA 
element  library  is  described  in  chapter  4. 

At  this  stage,  a  first  estimation  of  the  computation  size  can 

generally  be  made,  using  a  gross  determination  of  the  node  and 
degree  of  freedom  numbers  and  previous  CPU  time  measurements  or 
classical  numerical  rules.  This  estimation  is  important  since  It 
can  avoid  the  generation  of  unrealistic  meshes.  At  this  stage 
also,  the  user  has  to  gather  the  requested  physical  parameters,  as 
described  in  chapter  4,  and  to  ensure  their  accuracy,  which  very 
definitely  determines  the  accuracy  of  the  result.  Though  not 
directly  involved  in  the  data  file  generation,  this  stage  is 
fundamental . 

1.3.2.  THE  MESH  GENERATION. 

-Following  the  previous  choice  of  element  type(s),  as  defined 
during  the  first  stage,  the  mesh  generation  Implies  the  splitting 
of  the  whole  domain  into  a  given  number  of  elements.  This  Is 
realised  by  defining  nodes  and  assigning  them  to  corresponding 
elements.  Node  definition  is  made  by  listing  the  node  coordinates 
In  a  given  order,  called  the  node  numbering  order.  Then, 
assignment  of  the  nodes  to  the  corresponding  elements  is  obtained 
by  listing  for  each  element,  in  a  given  order  named  its  topology, 
the  nodes  which  belong  to  It.  The  node  coordinate  description  Is 
detailed  in  chapter  3,  together  with  all  other  entries.  The 
topology  description  is  presented  also  in  chapter  3,  after  the 
node  coordinate  description,  but  more  details  are  available  on  an 
element  by  element  basis  in  chapter  4. 
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During  the  mesh  generation,  the  user  can  use  the  isoparametric 
character  of  all  the  ATILA  elements,  which  allows  easy  curved  line 
modelling  and  lower  density  meshes.  Nevertheless,  a  great  deal  of 
attention  has  to  be  paid  to  the  element  geometrical  aspect,  as 
well  as  to  the  mesh  step  which  has  to  be  related,  for  example,  to 
the  minimal  wavelength  value  of  interest. 

Mesh  generation  can  be  fully  carried  out  by  the  user.  However, 
in  many  cases,  partial  regularity  of  the  mesh  allows  the  use  of  an 
automatic  generation  which  provides  node  coordinates  and  element 
topologies.  Then,  the  user  has  only  to  define  a  gross  splitting, 
built  with  super  elements,  and  to  select  their  automatic  splitting 
into  finite  elements.  This  easier  procedure,  described  in  chapter 
6,  relies  upon  the  use  of  the  MOSAIQUE  code. 

1.3.3.  THE  DATA  FILE  PREPARATION. 

An  ATILA  data  file  includes  entries,  specific  data  and 
comments.  It  is  devoted  to  the  description  of  the  type  of 
analysis,  the  node  coordinates,  the  element  topology,  the  material 
properties,  the  element  geometrical  properties,  which  are  all 
entries,  and  the  loading  and  boundary  conditions  which  are 
formatted  data.  In  some  cases,  it  also  includes  specific  mixed 
finite  element/plane  wave  model  data  and  graphic  display  data".  The 
whole  data  file  preparation  is  described  in  chapter  3  which  has  to 
be  carefully  read.  Specific  element  information  which  are  required 
during  this  preparation  can  also  be  found  in  chapter  4.  Finally, 
when  the  automatic  mesh  generation  code  MOSAIQUE  is  used,  data 
file  preparation  can  be,  in  a  large  part,  automatically  obtained, 
as  explained  in  chapter  6. 

A  detailed  check  of  the  data  file  can  be  processed  by  running  a 
specific  progam,  named  PGRAPH,  which  produces  a  complete  listing 
of  the  various  data,  Including  the  generated  degree  of  freedom 
numbering  which  reflects  exactly  the  prescribed  boundary 
conditions,  as  well  as  different  kinds  of  graphic  displays  (total 
or  partial  mesh,  with  or  without  element,  node  or  degree  of 
freedom  number Ings).  The  use  of  PGRAPH  is  described  in  chapter  5. 

1.3.4.  THE  RUNNING  OF  THE  CODE. 

Running  of  the  code  requires  the  edition  of  a  main  program 
which  has  to  set  the  array  sizes,  call  successive  subroutines  to 
read  the  data  file,  compute  the  elementary  stiffness  and.  mass 
matrices,  assemble  these  matrices  into  the  global  stiffness  and 
mass  matrices,  solve  the  equations  and  display  the  results.  This 
progam  generation  Is  automatic,  given  the  corresponding  data  file, 
and  is  provided  by  running  the  program  PGEN.  PGEN  will  issue 
successive  questions,  on  the  terminal  screen,  to  specify  various 
user's  requests,  such  as  array  storages  for  a  delayed  solving, 
array  storages  for  a  postprocessing,  graphic  displays...  It  will 
provide  various  files,  including  the  main  program  FORTRAN  source 
and,  in  the  case  of  a  VAX  computer  under  VMS  operating  system,  a 
COMMAND  file  to  compile  this  main  program  and  edit  the  links. 
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Running  of  the  code  will  provide  a  detailed  result  listing  and, 
according  to  the  specific  requests,  graphic  displays  or  stored 
arrays  for  restarting  another  solving  or  for  postprocessing.  It 
has  to  be  noted  that  PGEN  is  also  able  to  generate  specific  main 
programs  to  restart  given  runs.  Every  information  which  concerns 
the  program  running  as  well  as  PGEN  Is  contained  in  chapter  5. 

1.3.5.  RESULT  USING  AND  POSTPROCESSING. 

As  soon  as  the  run  Is  completed,  a  LISTING  file  Is  available 
which  contains  all  the  results:  eigenf requencies  and  eigenvectors, 
displacement  and  pressure  fields,  electrical  impedance...  If 
requested,  specific  arrays  are  also  stored  for  postprocessing,  and 
mainly  the  array  PLO  which  contains  displacement  and  pressure 
fields  for  every  loading  case,  mode  number  or  excitation 
frequency.  Graphic  displays  of  the  distorted  structure  can  be 
easily  obtained  on  various  workstations  using  PGRAPH .  -  PGRAPH  will 
issue  FORTRAN  and  COMMAND  files  which  rely  upon  the  GKS  package 
and  have  to  be  executed  to  provide  the  display.  As  mentioned 
previously,  the  use  of  PGRAPH  is  described  In  chapter  5. 

1.3. 6. SUMMARY. 

As  a  summary,  the  next  page  displays  a  simple  flow  chart  ^>f  an 
ATILA  finite  element  modelling.  A  large  part  of  this  flow  chart 
will  be  reproduced  in  a  more  detailed  and  accurate  fashion  in 
chapter  5,  with  reference  to  the  various  files  and  programs. 
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1.  MODEL  DEFINITION. 

Choice  of  the  type  of  analysis. 

Choice  of  the  elements. 

Choice  of  the  solving  algorithms. 

Computation  size  estimation. 

Gathering  of  all  physical  parameters. 


2.  MESH  GENERATION  (USER  OR  MOSAIQUE) 

Node  definition:  node  coordinate  list 

Node  appropriation  to  the  elements: 

element  topologies. 

Check  of  element  geometrical  aspects. 
Check  of  mesh  steps. 


3.  DATA  FILE  PREPARATION. 

Definition  of  entr I es, including  node 

coordinate  list  and  element  topologies. 

Definition  of  specific  data:  loads, 

boundary  conditions,  mixed  model. 

Check  of  the  mesh  and  data  file:  PGRAPH. 


4.  RUNNING  OF  THE  CODE  (FIRST  RUN  OR  RESTART) 

Generation  of  the  main  program:  PGEN. 

Specification  of  user's  requests:  array 

storages,  graphic  displays... 

Run  processing. 


5.  RESULT  USING  AND  POSTPROCESSING. 

Reading  of  the  LISTING  file. 
Graphic  display  of  the  distorted  shape  of 
the  structure:  PGRAPH. 
Postprocessing  using  various  resulting  arrays 
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CHAPTER  2 
TYPES  OF  ANALYSIS 


THE  FINITE  ELEMENT  CODE  AT  I  LA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  "INSTITUT  SUPERIEUR  D* ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D' ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN.  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE). 
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2.1  GENERAL  FORMULATION 


This  section  describes  the  general  formulation  which  has  been 
retained  to  model  radiating  piezoelectric  transducers  using  the 
ATILA  code  and  provides  a  complete  list  of  the  available  types  of 
analysis.  Then,  for  each  type  of  analysis,  the  following  sections 
give  more  information  and  specific  details. 


Boundary 
surface  S. 


To  model  a  radiating  piezoelectric  transducer  using  the  ATILA 
code,  the  finite  element  mesh  must  include  the  structure  as. well 
as  a  part  of  the  fluid  domain,  as  described  on  the  above  figure. 
Thus,  the  unknown  fields  are  the  displacement  field  u  in  the  whole 
structure,  the  electrical  potential  v  In  the  piezoelectric 
material  and  the  pressure  p  in  the  fluid.  Then,  the  equations 
which  are  solved  are,  first,  the  equation  of  motion  In  the 
structure,  second,  the  Poisson  equation  in  the  piezoelectric 
material  and,  third,  the  Helmholtz  equation  in  the  fluid. 
Kinematic  and  dynamic  continuity  equations  between  displacement 
and  pressure  fields  are  prescribed  at  the  interface,  due  to  the 
variational  formulation,  and  an  appropriate  damping  condition  is 
applied  to  the  spherical  external  surface  S  that  surrounds  the 
fluid  doma  i  n . 
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In  matrix  form,  the  complete  set  of  equations  associated  with 
this  problem  is: 


[Kuu]  -  w2[M] 

[Kuv  ] 

-a2c2w2[L]T 
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where: 

-  U,  F,  V,  Q  and  P  are  vectors  that  contain  the  nodal  values  of 
the  displacement  field,  the  lumped  applied  forces,  the 
electrical  potential  and  charges,  the  pressure  field, 

-  F*  is  a  vector  that  contains  the  nodal  values  of  the  outgoing 
flux  through  S, 

-  [Kuu]  and  [H]  are  the  solid  and  fluid  stiffness  matrices, 

-  [M]  and  [M-|]  are  the  solid  and  fluid  mass  matrices, 

-  [Kuv]  and  [Kvv]  result  from  the  assembling  of  the  piezoelectric 
and  dielectric  matrices, 

-  [L]  is  the  fluid-structure  interface  connectivity  matrix, 

-  [0]  is  a  nul I  matr Ix, 

-  a  and  c  are  the  fluid  density  and  sound  speed, 

-  w  is  the  circular  frequency, 

-  T  means  transposed. 

In  fact,  internal  electrical  potential  degrees  of  freedom  are 
associated  with  null  nodal  charges  and  can  be  easily  condensed. 
Moreover,  nodes  which  belong  to  hot  electrodes  have  the  same 
electrical  potential  degree  of  freedom  which  can  be  factored, 
following  simple  algebraic  operations.  Finally,  the  reference 
potential  can  be  arbitrarily  put  equal  to  zero.  Thus,  the 
preceeding  matrix  equation  Is  automatically  transformed  by  the 
ATILA  code  into  the  following  form: 


[Km,]  -  w2[M] 
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where  elastic,  piezoelectric  and  dielectric  terms  are  modified  by 
the  internal  electric  potential  static  condensation,  K  uv  and  0 
are  now  vectors,  Kvv  is  now  a  scalar,  V  is  the  applied  electrical 
potential  and  I  Is  the  electrical  current  entering  into  the 
structure.  Then,  if  a  dipolar  damping  condition  is  applied  on  the 
external  boundary  surface  S,  F'  is  given  by: 


F*  -  -(1/ac).(1/R  +  Jw/c).[D]  P 


+(1/ac).{(1/R  -  Jw/c)/(1  +  w2R2/c2)}.[D']  P 


Chapter  2.  Types  of  analysis. 


where  R  Is  the  radius  of  S,  [D]  and  [D']  are  obtained  by 
assembling  radiating  surface  finite  elements  from  S,  the  first 
term  is  the  monopolar  contribution  associated  with  a  spherical 
wave  impedance  and  the  second  term  is  the  dipolar  contribution. 

Starting  with  this  general  problem,  a  lot  of  different  specific 
analysis  can  be  processed.  Those  for  which  AT  I  LA  provides  a  solver 
are  the  fol lowing: 

-  computation  of  the  displacement  field  of  an  elastic  structure 
under  a  static  lumped  or  distributed  load  (STA  1), 

-  computation  of  the  displacement  field  and  the  electrical 
potential  for  a  piezoelectric  structure  under  a  static  lumped  or 
distributed  load  (STA  2), 

-modal   analysis   of   an   elastic   structure,   providing   the 

eigenf requencies  and  eigenmodes  (MOD  1), 
-modal  analysis  of  a  piezoelectric  structure,  providing  the 

eigenf requencies  and  eigenmodes  under  resonance  or  ant  I  resonance 

electrical  conditions  (MOD  2), 

-  modal  analysis  of  a  closed  fluid  domain,  with  zero  pressure 
and/or  zero  flux  boundary  conditions,  providing  the 
eigenf requencies  and  eigenmodes  (  MOD  3), 

-modal  analysis  of  a  closed  hydroelastic  system,  built  with  two 
coupled  elastic  and  closed  fluid  domains,  providing  the 
eigenf requencies  and  eigenmodes  (MOD  4), 

-  modal  analysis  of  an  elastic  structure  flooded  In  an 
incompressible  fluid  (added  mass  effect  only),  providing  the 
eigenf requencies  and  eigenmodes  (MOD  5), 

-  harmonic  analysis  of  an  in-air  piezoelectric  structure, 
providing  its  electrical  impedance  and  displacement  field  for  a 
given  excitation  frequency  (HAR  1), 

-  harmonic  analysis  of  an  in-air  piezoelectric  structure,  using 
the  mixed  plane  wave/finite  element  method  and  providing  the 
displacement  field  for  a  given  excitation  frequency  (HAR  2), 

-  radiated  pressure  computation  for  a  given  displacement  field  of 
a  vibrating  structure  at  a  given  frequency  (HAR  3), 

-harmonic  analysis  of  a  radiating  piezoelectric  transducer, 
providing  the  electrical  impedance  and  the  displacement  and 
pressure  fields  in  the  whole  mesh  for  each  given  excitation 
frequency  (HAR  4) , 

-  harmonic  analysis  of  a  radiating  piezoelectric  transducer,  using 
the  mixed  plane  wave/finite  element  method  and  providing  the 
displacement  and  pressure  fields  in  the  whole  mesh  for  each 
given  excitation  frequency  (HAR  5). 
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2.2  STATIC  ANALYSIS  OF  AN  ELASTIC  STRUCTURE  (STA  1) 


In  this  case,  only  the  displacement  field  is  concerned,  while  w 
has  to  be  put  equal  to  zero,  and  thus  the  whole  matrix  equation  Is 
reduced  to: 


CKUU]  U 


The  user  has  to  provide  the  loading  F  and  the  solver  computes  U.  F 
has  to  be  defined  strictly  for  the  finite  element  model,  which 
means  that,  for  example,  it  has  to  be  divided  by  2*pl  If  the  model 
Is  ax  isymmetr  leal  ,  by  2  if  it  Is  applied  in  a  symmetry  plane 
limiting  the  mesh...  A  frontal  assembling  Is  performed,  coupled 
with  a  static  Gauss  condensation,  and  then  a  classical 
backsubst I  tut  Ion  algorithm  restores  the  whole  response.  Several 
different  loading  cases  can  be  solved  simultaneously.  The  user  has 
to  take  care  of  the  presence  of  rigid  body  modes,  which  are 
associated  with  singularities  of  the  system  of  equations.  FiriaMy, 
if  distibuted  loads  are  applied,  the  user  has  to  compute  the 
corresponding  lumped  forces,  using  the  element  interpolation 
funct  ions. 
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2.3  STATIC  ANALYSIS  OF  A  PIEZOELECTRIC  STRUCTURE  (STA  2) 


In  this  case,  the  matrix  equation  is  reduced  to 


[Kuu]    K  UY 
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If  the  external  electrical  connection  is  a  short  circuit, 
equal  to  zero.  The  proceeding  equation  is  reduced  again  to: 


V  is 


[Kuu]  U  «  F 


The  user  has  to  provide  the  loading  F  and  the  solver  computes  U, 
following  the  same  procedure  as  in  the  proceeding  section.  If  the 
external  electrical  connection  is  open,  Q  Is  equal  to  zero  and: 


[Kuu]    !£  uv 
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kVV 
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The  user  has  to  provide  F  and  the  solver  computes  U  and  the 
induced  electrical  potential  V.  Here  again,  attention  has  to  be 
paid  to  rigid  body  modes  as  well  as,  in  the  second  case,  to  the 
uniqueness  of  the  elect  lea  I  potential.  As  for  STA  1,  several 
loading  cases  can  be  processed  at  the  same  time. 
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2.4  MODAL  ANALYSIS  OF  AN  ELASTIC  STRUCTURE  (MOD  1) 


In  this  case,  only  the  displacement  field  is  of  concern  and  the 
loading  force  is  put  equal  to  zero.  The  matrix  equation  is  reduced 
to: 


{[Km,]  -  w2[M]}  U 


Eigenvalues  and  eigenvectors  from  this  linear  system  of  equations 
are  the  eigenfrequencies  and  eigenmodes  of  the  structure. 

The  ATILA  code  uses  frontal  assembling  coupled  to  a  Guyan 
condensation.  Thus,  the  user  has  to  provide  the  number  of  master 
degrees  of  freedom  which  have  to  be  retained  in  the  final 
diagona I  I  sat  ion  phase.  To  do  this,  he  has  to  keep  in  mind  that  the 
CPU  time  which  is  reqested  for  the  diagona I i sat  ion  phase  Is 
roughly  proportional  to  the  cube  of  the  master  degree  number, 
absolute  value  of  this  time  depending  upon  the  algorithm,  the 
required  accuracy,  the  computer  and  the  compiler  option.  He  has 
also  to  keep  In  mind  that  accurate  eigenfrequencies  are  obtained 
only  for: 


w2  «  Min(K| | /M | |) 


where  K| |  and  MM  are  the  diagonal  terms  of  the  stiffness  and  mass 
matrices  associated  with  the  condensed  slave  degrees  of  freedom. 
Thus,  he  has  to  perform  a  compromise  between  these  opposing 
requirements.  For  a  given  master  degree  number,  the  choice  of  the 
specific  degrees  which  have  to  be  retained  can  be  provided  by  the 
user  or  automatically  selected  by  the  code,  using  the  Henschell's 
cr  iter  la. 

The  diagonal  I  sat  Ion  algorithm  which  is  used  by  ATILA  relies 
upon  the  classical  Householder  tr idiagona I isat ion  method,  coupled 
to  a  Sturm  sequence.  It  can  be  processed  in  simple  or  double 
precision.  CPU  times,  as  obtained  on  an  IBM  4341  Mod.  2  (4. MO), 
are  reported  on  the  graphic  displays  of  pages  8  and  9.  Recently,  a 
new  diagonal isat Ion  algorithm  has  been  implemented.  It  provides 
the  eigenfrequencies  in  double  precision  and  the  eigenvectors  in 
single  precision.  Though  much  more  efficient,  it  is  still  under 
test  and  can  be  only  used  by  advanced  users,  following  a  short 
modification  of  the  ATILA  main  program*.  In  all  the  cases,  ATILA 
provides  all  the  eigenfrequencies,  the  number  of  which  is  equal  to 
the  number  of  master  degrees,  and  only  the  first  twenty 
eigenmodes,  if  available.  It  provides  also  the  rigid  body  modes, 
which  belong  to  the  previously  mentioned  first  twenty  modes.  Upon 
the  user's  request,  a  backsubst I  tut  Ion  routine  can  provide  the 
slave  degrees  of  freedom  for  any  of  the  obtained  eigenmodes. 
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In  some  cases,  care  has  to  be  taken  because  of  mixing  between 
modes  that  have  frequencies  which  are  too  close,  as  well  as  of 
matrix  i I l-cond 1 1 lonlng.  Successive  diagonal  I  sat  Ion  with  different 
precision  or  different  master  degree  numbers  can  illuminate  these 
two  points. 


*.  To  use  the  new  diagonal  I  sat  ion  algorithm,  the  user  has  to 
edit  the  main  program  FORTRAN  source,  as  obtained  for  a  double 
precision  eigenvalue  computation,  replace  CALL  D52200  (DLIST,S, 
S,CPDDC,PLO, IP,DM,Z)  by  CALL  S52200  (DL IST,S,CPDDC,PLO, IP  ,Z)  and 
delete  the  DM  array  from  the  DIMENSION  set  of  instructions. 
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CPU  TIME  FOR  THE  ASSEMBLING  PROCESS  WITH  RESPECT  TO  THE  NUMBER  OF 
MASTER  DEGREES  OF  FREEDOM.  The  parameter  Is  the  total  number  of 

degrees  of  freedom. 
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CPU  TIME  FOR  THE  D I AGONAL  I  SAT  I ON  PROCESS  WITH  RESPECT  TO  THE 
NUMBER  OF  MASTER  DEGREES  OF  FREEDOM.  The  parameter  is  the  total 

number  of  degrees  of  freedom. 
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2.5  MODAL  ANALYSIS  OF  A  PIEZOELECTRIC  STRUCTURE  (MOD  2) 


In  this  case,  the  matrix  equation  is  reduced  to: 


[Km.]  -  w2[M] 


£  uv 


K  uv 

1      1     ■ 

!      0 

Kvv 

1   »  , 

1-l/jw 

f  a  resonance  condition  occurs,  V  is  equal  to  zero  and 


C[KUQ]  -  w2[M]}  U  -  0 


If  an  ant i resonance  condition  occurs,  I  Is  equal  to  zero  and: 


{[Kuu]  -  w2[M]  -  K  UV.K  UVT/KVV}  U  -  0 


In  the  two  cases,  all  the  information  given  in  the  proceeding 
section  with  respect  to  the  modal  analysis  of  an  elastic  structure 
is  applicable.  The  only  point  which  has  to  be  added  is  that,  due 
to  the  automatic  static  condensation  of  the  electrical  potential 
degrees  of  freedom,  all  these  degrees  of  freedom  are  slaves. 


.10. 


Chapter  2.  Types  of  analysis 


2.6  MODAL  ANALYSIS  OF  A  CLOSED  FLUID  DOMAIN  (MOD  3) 


In  this  case,  only  the  pressure  field  Is  concerned.  Moreover, 
as  zero  pressure  and/or  zero  flux  boundary  conditions  have  to  be 
applied,  the  F*  components  which  are  able  to  appear  in  the 
equation  are  all  zero.  Thus: 


in  section  2.5  with  respect  to  the 
structure  Is  applicable.  Moreover, 


{[H]  _  wz[M.,j)  P  -  0 

All  the  information  provided 
modal  analysis  of  an  elastic 
the  user  has  to  remember  that  the  zero  flux  condition  is  the 
natural  boundary  condition  of  this  finite  element  formulation  and 
automatically  applies  to  any  boundary,  if  not  superseded  by 
another  one.  This  modal  analysis  can  be  used  to  determine  the 
irregular  frequencies  associated  with  an  exterior  Helmholtz 
integral  equation. 


1  1 
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2.7  MODAL  ANALYSIS  OF  AN  HYDROELASTIC  SYSTEM  (MOD  4) 


In  this  case,  structure  and  fluid  domains  exist  together  and 
interact  along  their  interface.  However,  the  fluid  domain  has  to 
be  closed,  no  radiation  damping  being  included  In  the  model.  Thus, 
fluid  surfaces  which  are  not  in  contact  with  the  solid  structure 
have  to  be  zero  pressure  or  zero  flux  surfaces,  as  Is  the  case, 
for  example,  for  fluid  in  an  open  tank.  Restricting  the  general 
matrix  equation  of  section  2.1  to  the  displacement  and  pressure 
field,  the  equation  which  Is  obtained  is: 


[Kuu]  -  w2[M]     -  [L] 
-aac2w2[L]T   [H]  -  w2[M-|] 


:u: 


i  p  i 


Apart  for  a  slight  difference  regarding  the  algorithm,  all  the 
information  given  in  section  2.4  with  respect  to  the  modal 
analysis  of  an  elastic  structure  is  applicable.  With  respect  to 
the  algorithm,  the  above  matrix  has  to  be  made  symmetrical  before 
the  AT  I  LA  classical  diagona I isat ion  is"  performed.  This 
symmetr Isat Ion  is  processed  using  a  standard  method  due  to  Irons. 
However,  it  has  to  be  pointed  out  that  the  symmetr  Isat  Ion  method 
requires  the  inversion  of  [H]  and  then  implies  that  no  fluid 
isopressure  mode,  equivalent  to  structure  rigid  body  modes, 
exists.  Thus,  at  least  a  small  part  of  the  fluid  surface  has  to  be 
a  zero  pressure  surface.  Recent  developments  have  led  to  the 
design  of  a  new  solver  which  avoids  this  difficulty,  but  It  is 
still  under  test  and  can  only  be  run  by  advanced  users  .  Moreover, 
during  the  assembling  phase,  as  well  as  for  an  elastic  structure, 
Guyan  condensation  is  performed  for  both  solid  and  fluid  degrees 
of  freedom.  Care  has  to  be  given  to  the  interface  degrees  of 
freedom  which  have  to  be  kept  in  sufficient  number.  Finally, 
piezoelectric  parts  can  be  included  in  the  structure.  Then,  the 
user  has  also  to  provide  the  corresponding  electrical  boundary 
condi t  ions. 


*.  To  use  this  new  algorithm,  the  user  has  to  edit  the  main 
program  FORTRAN  source,  as  obtained  for  this  analysis,  delete  CALL 
R52300  and  replace  CALL  R52200  by  CALL  S52300  with  the  same 
parameters.  The  new  algorithm  detects  automatically  the 
isopressure  or  rigid  body  modes  and  deals  with  them  using  a 
deflation  method.  In  opposition  to  the  standard  algorithm,  it  can 
also  compute  the  eigenmodes  for  mere  fluid  or  solid  domains, 
without  interaction. 
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2.8  MODAL  ANALYSIS  OF  AN  ELASTIC  STRUCTURE  FLOODED  IN  AN 
INCOMPRESSIBLE  FLUID  (MOD  5) 


Incompressible  fluid  condition  is  an  approximation  that  is 
useful  to  compute  the  fluid  mass  loading  effect  on  the  low 
frequency  modes  of  an  elastic  structure.  Obviously,  It  discards 
any  radiation  effect.  If  It  is  applied  to  the  matrix  equation  of 
the  previous  section,  this  equation  becomes: 

{[Kuu]  -  w2([M]  +  a2c2[L][H]-1[L]T)}  U  -  0 

where  a2c2[L][H]~1 [L]T  is  the  added  mass.  Then,  'the  classical 
ATILA  diagonal i sat  ion  provides  the  structure  eigenf requencles  and 
eigenmodes.  To  obtain  this  particular  condition  using  ATILA 
requires  that  all  the  fluid  degrees  of  freedom  are  condensed,  this 
approximation  producing  exactly  the  same  result  as  the 
incompress ibi I i ty  hypothesis  .  If  the  fluid  domain  extension  is 
infinite,  it  has  to  be  simply  truncated  before  the  mesh 
generation.  Successive  truncations  at  various  distances  from  the 
wet  structure  are  the  best  convergence  test. 


*.  After  the  Guyan  condensation  of  all  the  fluid  degrees  of 
freedom,  the  matrix  which  has  to  be  diagonal ised  is  already 
symmetrical  but  the  automatic  program  generator  will  provide  the 
call  for  the  symmetr isat ion  subroutine  because  the  data  file 
corresponds  to  a  fluid-structure  modal  analysis.  Thus,  before 
running  ATILA  the  user  must  edit  the  FORTRAN  file  and  delete  the 
I ine  CALL  R52300. 
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2.9  HARMONIC  ANALYSIS  OF  AN  IN-AIR  PIEZOELECTRIC  STRUCTURE  (HAR  1) 


In  this  case,  an  electrical  potential  of  given  frequency  Is 
applied  to  an  in-air  free  piezoelectric  structure.  Then,  the 
general  matrix  equation  of  section  2.1  can  be  reduced  to: 


[K,1U]  -  w2[M] 


kUU 


■i  uv 


!S  uv 


"VV 


:u: 


!V! 


0 

-l/Jw 


To  obtain  this  whole  equation  using  the  ATILA  code  implies  that 
the  applied  electrical  potential  degree  of  freedom  is  not 
automatically  condensed,  which  requires  a  specific  identification 
using  the  boundary  condition  data.  Then,  the  solver  computes  the 
displacement  field  U,  using  the  first  line  of  the  above  matrix 
equat  ion: 


{[Kuu]  - 


w2[M]} 


-  -K 


uv 


.V 


as  well  as  the  electrical  Impedance  Z,  using  the  second  line 


1/Z  -  -jw  (Kvv  +  K  UVT.U/V) 


The  solving  algorithm  is  a  simple 
run,  solving  can  be  performed 
frequencies.  Moreover,  job  restart 
values  without   performing   again 
impedance  value  which  is  provided 
finite  element  model.  Thus,  the  user 
2*pi  if  the  model  Is  ax  Isymmetr  leal 


Gauss  elimination.  In  the  same 
for  20  different  excitation 

is  possible  for  new  frequency 

the  whole   assembling.   The 

corresponds  strictly  to  the 

has  to  divided  the  result  by 

2  if  the  mesh  is  restricted 


to  half  of  the  structure  by  a  symmetry  plane,  etc 
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2.10  HARMONIC  ANALYSIS  OF  AN  IN-AIR  PIEZOELECTRIC  STRUCTURE  USING 

THE  MIXED  METHOD  (HAR  2) 


In  this  case,  the  driving  part  of  the  structure  is  modelled 
using  the  plane  wave  approximation  and  a  transfer  matrix  method, 
while  the  remaining  Inactive  part  of  this  structure  Is  modelled 
using  finite  element.  These  two  parts  are  in  contact  along  a  plane 
interface,  which  exhibits  also  the  plane  wave  type  displacement 
field.  Thus,  all  the  degrees  of  freedom  which  belong  to  the  finite 
element  side  of  the  interface  and  correspond  to  a  mechanical 
displacement  normal  to  it  have  to  be  identical  and  are  denoted 
hereafter  as  U0.  Then,  the  transfer  matrix  method  provides,  for 
each  given  frequency,  a  mechanical  force  F0  which  is  generated  by 
the  piezoelectric  effect  and  is  applied  to  the  finite  element  part 
at  the  interface.  This  force  can  be  written: 


-  Z0.U0  +  Fv 


where  Z0  is  a  mechanical  impedance  and  Fv  contains  the  electrica 
excitation  effect.  Finally,  the  general  matrix  equation  becomes: 


[KQU]  -  w2[M]   K 


£  O 


!i  O 

Koo 


U 


where  [Kuu]  and  [M]  correspond  to  all  the  degrees  of  freedom  apart 
from  U0.  The  user  has  to  select  the  excitation  frequencies  and  to 
specify,  using  the  correct  boundary  condition  data,  the  degree  of 
freedom  U0.  The  solving  algorithm  is  a  simple  Gauss  elimination, 
and  the  code  provides  the  displacement  field  in  the  finite  element 
part  of  the  structure.  Twenty  different  excitation  frequencies  can 
be  provided  for  the  same  run.  Moreover,  as  in  section  2.9,  job 
restart  is  possible  with  new  frequency  values.  This  mixed  method 
is  generally  well  suited  to  the  modelling  of  Tonpilz  type 
transducers  and  allows  large  savings  of  CPU  time  . 
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2.11  RADIATED  PRESSURE  COMPUTATION  FOR  A  GIVEN  DISPLACEMENT  FIELD 

(HAR  3) 


In  this  case,  a  vibrating  structure  radiates  an  acoustic  wave 
in  an  infinite  medium  and  the  vibration  displacement  is  assumed  to 
be  known  for  each  given  frequency.  The  general  matrix  equation 
from  section  2.1  can  be  reduced  to: 


£[H]  -  w2[M-|]}  P  -  ac2F'  +  a2c2w2[L]T.U 


Then,  if  monopolar  or  dipolar  damping  is  used,  following  the  third 
equation  of  section  2.1,  the  damping  term  can  be  written  under  the 
form: 


ac^F*  -  [G]  P 


where  [G]  stands  for  a   linear  complex  frequency  dependent 
operator.  Thus: 


{[H]  -  w2[M-,]  -  [G]}  P  -  a2c2w2[L]T.U 


The  user  has  to  define  U  and  the  given  excitation  frequencies. 
To  define  U,  he  has  to  write  a  small  subroutine,  named  FABU* , 
which  contains  the  displacement  value  for  each  interface  solid 
degree  of  freedom,  and  has  to  be  added  to  the  main  program.  Only 
one  displacement  field  can  be  defined  at  a  time  and  this  field  has 
to  be  real  (phase  reference).  Ten  successive  frequencies  can  be 
processed  in  the  same  run  and,  as  for  the  previous  section, 
restart  is  available.  The  solver  splits  the  previous  complex 
equation  into  two  real  equations.  Then,  a  matrix  Inversion 
followed  by  a  Gauss  elimination  is  performed.  For  each  frequency, 
ATILA  provides  the  real  and  imaginary  parts  of  the  pressure  field, 
printed  in  two  successive  tables. 

To  obtain  accurate  results,  the  user  has  to  give  a  great  deal 
of  attention  to  the  damping  element  choice.  If  monopolar  damping 
is  used,  the  external  boundary  surface  has  to  be  centered  at  the 
acoustic  center  and  it  has  to  be  in  the  farfield.  If  dipolar 
damping  is  used,  the  external  boundary  surface  has  also  to  be 
centered  at  the  acoustic  center  but,  generally,  it  can  be  deeply 
in  the  nearfield.  However,  in  this  second  case,  farfield 
information  is  only  available  after  an  appropriate  postprocessing. 
Dipolar  elements  can  only  be  used  for  ax isymmetr ica I  models. 
Three-dimensional  dipolar  element  design  Is  In  progress. 
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Degrees  of  freedom  which  belong  to  the  exterior  boundary 
surface  are  automatically  master  degrees.  Those  which  belong  to 
the  interface  can  be  slave  degrees  but  the  user  has  to  keep  a 
sufficient  number  of  them  as  master  degrees.  Finally,  care  has  to 
be  taken  to  keep  the  maximal  distance  between  two  successive  fluid 
master  nodes  in  the  propagation  direction  less  than  one  quarter  of 
a  wavelength. 


*.  The  FORTRAN  subroutine  which  has  to  be  added  at  the  end  of 
the  main  program,  to  provide  the  given  displacement  field  U,  is 
FABU  (CPDDC,U, IOLE, IP).  It  is  called  by  the  solving  subroutine 
R52400  during  the  run.  IOLE  and  IP  are  control  Integers  which  are 
defined  In  R52400,  as  well  as  the  CPDDC  array.  In  return,  FABU 
provides  all  the  components  of  U,  numbered  in  the  ascending  order. 
U  has  to  contain  IDT  components,  where  IDT  is  the  total  number  of 
degrees  of  freedom,  but,  in  fact,  only  the  mechanical  degrees  of 
freedom  are  concerned.  The  following  lines  are  provided  .as  a 
standard  example: 

SUBROUTINE  FABU(CPDDC,U, IOLE , IP) 
DIMENSION  CPDDCUP,*),  U(*) 
U(1)  -  ... 
U(2)  -  ... 


RETURN 
END 
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2.12  HARMONIC  ANALYSIS  OF  A  PIEZOELECTRIC  RADIATING  STRUCTURE 

(HAR  4) 


This  is  the  general  case.  Using  the  operator  [G]  defined  In  the 
previous  section  2.11,  the  matrix  equation  is: 


[Km.]  ~  w2[M] 


£  uv 


-a2c2w2[L]T 


K.  uv 


kVV 

0 


-  [L] 


[H]  -  w2[M-|]|  IP 


:  :y 

!  0   ! 

!  !V 

-  l-l/Jwj 

!  !P 

!  [G]P| 

Then,  it  can  be  transform  into: 


[Kuu]  -  w2[M]         -  [L] 
-a2c2w2[L]T   CCH]  -  w2[M-,]  -  [G]) 


:  :y: 

i  i  i 
i  i  i 

_ 

-K  UV.V 

:  !P! 

0 

and: 


1/Z 


-jw 


_  uv  *y  "^  vv ' 


V) 


The  user  provides  the  excitation  frequencies  and,  using 
specific  boundary  condition  data,  defines  the  degree  of  freedom 
that  is  associated  with  the  applied  electrical  potential.  The 
solver  computes,  using  matrix  inversion  and  Gauss  condensation, 
the  real  and  imaginary  parts  of  the  displacement  and  pressure 
fields.  Then,  it  computes  the  electrical  impedance.  Real  and 
imaginary  parts  are  provided  in  two  successive  tables.  Ten 
successive  excitation  frequencies  can  be  processed  in  the  same 
run,  and,  as  in  the  previous  section,  restart  is  available. 

Information  given  in  the  previous  section  and  regarding  damping 
elements,  nearfleld  or  farfield  conditions,  boundary  surface  and 
interface  degrees  of  freedom  remain  available  for  this  section  and 
are  not  duplicated.  Information  given  in  section  2.9  regarding  the 
electrical  impedance  is  also  valid. 
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2.13  HARMONIC  ANALYSIS  OF  A  PIEZOELECTRIC  RADIATING  STRUCTURE 
USING  THE  MIXED  METHOD  (HAR  5) 


This  case  Is  a  mixing  of  the  previous  cases  reported  in 
sections  2.10  and  2.12.  Using  the  corresponding  notations,  the 
matrix  equation  becomes: 


[Km.]  -  w2[M] 


WUU 


£  O 


-a2c2w2[L]T 


Ko 


Koc 


-  CL] 


[H]  -  w2^] 


y 

:  o    : 

Uo 

-    |    FQ    | 

P 

!CG]P| 

The  user  provides  the  excitation  frequencies  and,  using 
specific  boundary  condition  data,  defines  the  degree  of  freedom 
that  is  associated  with  the  normal  displacement  to  the  interface. 
The  solver  computes,  using  matrix  inversion  and  Gauss 
condensation,  the  real  and  imaginary  parts  of  the  displacement  and 
pressure  fields.  Real  and  imaginary  parts  are  provided  in  two 
successive  tables.  Ten  successive  excitation  frequencies  can  be 
processed  In  the  same  run,  and,  as  in  the  previous  section, 
restart  is  aval l able. 

Information  given  In  section  2.11  regarding  damping  elements, 
nearfield  or  farfield  conditions,  boundary  surface  and  interface 
degrees  of  freedom  remain  available  for  this  section  and  are  not 
dupl icated. 
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CHAPTER  3 
DESCRIPTION  OF  A  DATA  FILE 


THE  FINITE  ELEMENT  CODE  AT  I  LA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  " INSTITUT  SUPERIEUR  D' ELECTRON  I QUE  DU  NORD" 
(I  SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D* ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN.  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE). 
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Chapter  3.  Description  of  a  data  file 


3.1  INTRODUCTION 


This  section  is  placed  first  in  the  data  preparation  part  of 
this  manual  for  pedagogical  purposes  only.  In  fact,  In  most  cases, 
an  automatic  data  generation  can  be  processed  using  the  MOSAIQUE 
code  and  a  more  simple  initial  data  file.  This  automatic  data 
generation  is  described  in  the  chapter  6.  When  used  first, 
MOSAIQUE  provides  data  files  which  are  identical  to  the  files 
described  hereafter. 

In  the  general  case,  an  ATILA  data  file  includes  five  different 
sect  ions: 

-  the  header,  which  contains  only  comments  and  describes  the  job 
under  consideration, 

-  a  set  of  entries,  which  allows,  among  other  things,  to  describe 
the  mesh  and  the  properties  of  the  materials, 

-  a  set  of  data,  which  describes  the  loading  and  boundary 
condi t  ions, 

-  a  set  of  data,  which  describes  the  plane  wave  part  of  the 
structure  In  the  case  of  a  mixed  model, 

-  a  set  of  data,  which  defines  the  graphic  displays. 

All  data,  except  for  the  angles,  must  be  given  using  MKS  units  or 
any  coherent  unit  system  which  is  deduced  from  the  MKS  system.  The 
angle  values  are  always  given  in  degrees. 

In  this  chapter,  we  always  refer  to  a  simple  illustrative 
example  and  give  the  corresponding  version  of  each  entry  or  data 
set.  This  example  is  the  modal  analysis  of  an  elementary 
ax isymmetr ical  piezoelectric  structure,  made  up  of  two  ceramic 
rings  and  two  thick  electrodes,  under  electrical  ant iresonance 
conditions.  This  structure  has  a  symmetry  plane  which  is  normal  to 
the  axis  and  the  modes  which  have  to  be  computed  are  symmetrical 
wijth  respect  to  this  plane.  Five  full  eigenvectors  are  requested. 
The  structure,  its  material  properties  and  the  corresponding  mesh 
are  described  on  the  following  two  pages.  Other  examples  are  given 
at  the  end  of  the  manual. 
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Test  example  description, 
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Dimensions  are  in  centimeters 


Material  properties. 


AU4G  ( a  I  urn  i  n I  urn  alloy), 

.Young's  modulus:  0.714E+11  Pa 
.Poisson's  ratio:  0.344 
.Density:       2780  kg/m3 


X5105  (piezoelectric  ceramic  material). 


Polarization:  axial  (up  for  the  top  ceramic  ring  and 

down  for  the  bottom  ceramic  ring). 
Elastic  constants  (MKS  units,  constant  electric  field) 

S11:   0.114E-10    S33:   0.126E-10 

S12:  -0.339E-11     S13:  -0.410E-11 

S44:   0.510E-10 
Piezoelectric  constants  (MKS  units): 

d33:   0.208E-09    d31 :  -0.955E-10 

d15:   0.491E-09 
Dielectric  constants  (MKS  units,  constant  strain): 

EP11:  0.667E-08    EP33:  0.687E-08 
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Mesh. 
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Due  to  simple  symmetry  arguments,  the  mesh  is 
of  the  structure.  Thus,  the  0Y  axis  is  the 
symmetry  plane  with  the  plane  of  the  figure. 


I imi ted  to  one  half 
intersection  of  the 
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3.2  THE  HEADER 


The  header  includes  comments  which  describe  the  Job  under 
consideration.  The  whole  header  is  automatically  copied  on  the 
first  page  of  every  listing  which  Is  related  to  the  job.  The  first 
line  of  the  header,  generally  a  title,  Is  automatically  repeated 
ahead  of  each  listing  page.  To  define  the  header,  each  line  begins 
with  the  symbol  *.  The  number  of  lines  in  a  header  is  unlimited. 


EXAMPLE 


+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
+  *...  MODAL  ANALYSIS  .  PIEZOELECTRIC  STRUCTURE  . 

+  * .  .  . 

+  * . . .  ANT  I  RESONANCE  CASE  

+  * . . .  AX  I SYMMETR I CAL  MODEL  

+  *...  1/2  OF  THE  STRUCTURE  (SYMMETRY  PLANE)  

+  * 

+  * . . .  NUMBER  OF  ELEMENTS :  2  

+  *...  NUMBER  OF  NODES:  13  

+  * . . .  NUMBER  OF  DOF :  26  

+  * . . .  NUMBER  OF  MASTER  DOF :  23  

+  *. . .  DOUBLE  PRECISION  

+  * 

+  * . . .  DATA  FILE  WR I TTEN  BY : 

+  * . . .  DATE  

+ 
+ 
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3.3  THE  ENTRIES 


The  various  entries  which  are  described  hereafter  are  listed  in 
a  given  order,  which  corresponds  to  a  classical  use  but  is  not 
mandatory.  These  entries  have  to  be  written  using  a  simple  free 
format,  described  below,  which  allows  an  easy  reading  of  the  file. 
An  entry  can  be  associated  to  given  peculiar  words,  named 
attributes,  which  have  to  be  written  on  the  entry  line,  or  to 
given  numerical  values,  named  parameters,  which  have  to  be  written 
on  the  following  line.  Entries  being  read  in  free  format,  all  the 
attributes  and  parameters  have  to  be  provided,  even  if  they  are 
dummy  or  if  their  numerical  value  is  zero. 


3.3.1  DESCRIPTION  OF  THE  FREE  FORMAT. 

Control  characters. 

.  [  ]  and  [,]  separate  two  words  on  the  same  line.  They  have  to  be 

used  to  separate  an  entry  from  its  following  attribute,  to 

separate  two  following  attributes  or  to  separate  two  following 

parameters. 
.  [/]  and  [-]  separate  two  following  lines  (carriage  return).  As 

well  as  the  carriage  return,  they  have  to  be  used  to  separate  an 

entry  from  Its  following  parameter. 
.  [?]  deletes  the  preceeding  characters  of  the  line, 
.  [*]  when  typed  at  the  beginning  of  a  new  line,  means  that  this 

line  is  a  comment.  This  line  is  not  printed  on  a  listing.  Such  a 

comment  has  not  to  be  confused  with  the  header. 
.  [* *]  allows  the  writing  of  a  comment  In  an  entry 

line.  This  comment  is  not  printed  on  a  listing. 

Real  number  syntax. 

The  following  syntaxes,  given  as  examples,  are  available: 

1.22E+02  .    -1.22E2 

-1E-3  .    1.22 

.  ~"  -1.  -1 

Examples. 

The  following  sequences  are  equivalent  to  define  the  given 

example:  one  alphanumeric  datum  (alpha)  on  the  first  line,  two 

integers  (1,2)  on  the  second  line,  three  real  numbers  (1.,2.,3.) 
on  the  third  I ine: 

.  alpha         .  alpha-1,2         .  alpha/1  2-1., 2.  3. 
1  ,2  1.2. ,3. 

1.,2.,3. 
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3.3.2  THE  ENTRY  DIRECTORY. 


NAME 

COMP. 

ATTRIBUTES  OR   PARAMETERS 

P 

ANAL YS 1 S 

YES 

STATIC,  MODAL  or  HARMONIC 

7 

CLASS 

NO 

AX  1  SYMMETRICAL,  PSTRESS  or 
PSTRAIN 

7 

NLOAD 

YES 

NLO 

7 

PRECISION 

NO 

SINGLE  or  DOUBLE 

8 

LCPDDC 

YES 

NLC 

8 

NODES 

YES 

node  coordinates 

(end  of  the  set:  dummy  line) 

9 

NEWAXES 

NO 

CARTESIAN,  CYLINDRICAL  or 

SPHERICAL 

XO  YO  ZO  OZZ  OYY  OXX 

10 

SCALE 

NO 

EX  EY  EZ 

11 

ELEMENTS 

YES 

Rxxxxx  MATERx  NGx 

Element  topology 

(end  of  the  set:  dummy  line) 

12 

MATERIALS 

YES 

MATERx 

Properties  of  the  material. 

(end  of  the  set:  dummy  line) 

14 

GEOMETRY 

NO 

NGx 

Geometrical  properties. 

(end  of  the  set:  dummy  line) 

17 

GEOMETRY 
POLAR IZAT. 

NO 

CARTESIAN,  CYLINDRICAL  or 

SPHERICAL. 

NGx 

Euler  angles,  coordinates. 

(end  of  the  set:  dummy  line) 

18 

REDUCTIONS 

YES 

1  N  MDF 

24 

FREQUENCIES 

NO 

Frequency  values. 

25 

RADIATION 

NO 

MONOPOLAR  or  DIPOLAR. 

25 

PRINTING 

NO 

NPR 

25 

END 

YES 

26 

COMP.  means  that  the  entry  Is  compulsory. 
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3.3.3  DESCRIPTION  OF  THE  ENTRIES. 


In  this  section,  the  use  of  the  symbol  6  In  front  of  the  entry 
name  means  that  this  entry  is  compulsory  and  has  to  be  used  in 
every  job.  The  other  entries  have  to  be  used  only  in  specific  jobs 
or  are  opt  ional . 


6   ANALYS I S  STATIC  or  MODAL  or  HARMONIC 

This  entry  defines  the  kind  of  analysis  which  is  requested  in  the 
job  (see  chapter  2,  theoretical  background)  and  is  compulsory. 


Test  example.     + 

+  ANALYSIS   MODAL 

+ 


CLASS  AX  I  SYMMETRICAL  Or  PLSTRESS  or  PLSTRAIN 

This  entry  specifies,  If  necessary,  that  the  analysis  corresponds 
to  an  ax isymmetr ical ,  a  plane  stress  or  a  plane  strain  model.  It 
is  only  needed  If  2D  elastic  elements  are  used,  but  then  is 
compulsory. 


Test  example.     + 

+  CLASS  AX  I SYMMETR I CAL 

+ 


6   NLOAD 
NLO 

This  entry  is  compulsory  and  provides  NLO  which  is: 

.  the  number  of  loading  cases  for  a  static  analysis, 

.  ^the  number  of  full  eigenvectors  which  are  required  In  the- case 

of  a  modal  analysis, 

.  the  number  of  given  frequency  values  in  the  case  of  an  in-air 

harmonic  analysis, 

.  twice  the  number  of  given  frequency  values  in  the  case  of  an 

harmonic  analysis  of  a  radiating  structure. 


Test  example.     + 

+  NLOAD 

+ 
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PRECISION  SINGLE  or  DOUBLE 

This  entry  is  only  used  if  an  eigenvalue  computation  is  run  after 
the  assembling  and  specifies  the  precision  of  this  computation. 
The  default  answer  is  single. 


Test  example.     + 

+  PRECISION  DOUBLE 

+ 


6   LCPDDC 
NLC 


The  CPDDC  array  contains  the  node  coordinates  (X,  Y,  Z)  and  the 
numbers  of  the  node  degrees  of  freedom  (translations  UX,  UY,  UZ, 
rotations  NX,  NY,  NZ,  pressure  P  and  electrical  potential  V) . 
These  data  are  listed,  with  one  row  per  node,  In  the  following 
order : 

X,  Y,  Z,  UX  (or  P),  UY,  UZ,  NX  (or  V),  NY,  NZ 

The  LCPDDC  entry  is  compulsory  and  provides  NLC,  which  Is  the 
number  of  terms  needed  In  one  row.  Thus,  in  the  case  of  an 
ax isymmetr ica I  analysis  of  an  elastic  structure  NLC  equals  5, 
while  in  an  ax isymmetr leal  analysis  of  a  piezoelectric  structure 
NLC  equals  7  and  in  an  analysis  of  a  plate  NLC  equals  8.  If 
useless  intermediate  columns  are  introduced,  due  to  the  given 
value  of  NLC,  the  corresponding  degrees  of  freedom  will  have  to 
be  deleted  later.  This  case  happens,  as  an  example,  for  the 
degrees  of  freedom  Included  in  columns  3,  4  and  5  in  the  plate 
case. 


Test  example.     + 

+  LCPDDC 
+  7  

+ 
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« 


NODES 
X1  Y1  Z1 
X2  Y2  Z2 
X3  Y3  Z3 


XN  YN  ZN 
.DUMMY.. 

This  entry  provides  the  node  coordinates,  listed  with  one  row  per 
node,  in  the  node  number  ascending  order.  The  list  has  to  be  ended 
by  a  dummy  line  and  the  entry  Is  compulsory.  If,  for  any  reason, 
the  node  coordinate  list  is  split  in  several  successive  sets,  each 
set  requires  a  new  entry  NODES  and  has  to  be  finished  by  a  new 
dummy  I Ine. 


WARNING 


When  the  model  includes  both  structure 
and  fluid  domains,  nodes  have  to  be  set 
In  the  following  order: 

SOLID  -  FLUID  -  RADIATING. 

For  ax  I  symmetrical  models,  the  OX  axis 
Is  necessarily  the  symmetry  axis. 


Test  example. 


+ 

+   NODES  

+0.00  0.02 

+0.00  0.05 

+0.00  0.08 

+  0.01  0.02 

+0.01  0.08 

+  0.02  0.02 

+   0.02  0.05 

+0.02  0.08 

+0.04  0.02 

+0.04  0.08 

+0.06  0.02 

+  0.06  0.05 

+  0.06  0.08 

+    

+ 
+ 
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NEWAXES  CARTESIAN  or  CYLINDRICAL  or  SPHERICAL 
XO  YO  ZO  OZZ  OYY  OXX 

The  Initial  coordinate  system  in  which  the  node  coordinates  are 
given  (default  system)  Is  cartesian  and  its  origin  and  axes  are 
the  absolute  origin  and  axes.  The  entry  NEWAXES  allows  it  to  be 
mod  I f  i  ed . 

If  the  new  system  Is  cylindrical  (spherical),  the  node  coordinates 
X,  Y,  Z,  associated  to  the  following  NODES  entry(ies)  become  the 
cylindrical  (spherical)  coordinates  R,  THETA,  Z  (R,  THETA,  PHI). 
Moreover,  XO,  YO,  ZO  are  the  coordinates  of  the  origin  of  the  new 
system,  while  OZZ,  OYY,  OXX  are  the  angle  values  of  the  Euler 
rotations  which  transform  the  initial  axes  In  the  new  axes.  The 
first  rotation  OZZ  is  around  the  OZ  axis  of  the  initial  system. 
The  second  rotation  OYY  is  around  the  OY  axis  provided  by  the 
first  rotation.  The  third  rotation  OXX  is  around  the  OX  axis 
provided  by  the  second  rotation.  Angle  values  have  to  be  given  in 
degrees. 

This  entry  can  be  used  several  times  during  the  description  of  the 
mesh,  its  parameters  being  always  defined  with  respect  to  the 
absolute  system.  To  come  back  to  the  absolute  system,  it  has  to  be 
used  with  the  CARTESIAN  attribute,  the  parameters  being  equal  to 
zero. 


Test  example. 


+ 

+   NODES    

+0.00      0.02    

+0.00     0.05    

+  0 .  00     0 .  08    

+  0.01      0.02    

+0.01      0.08    

+0.02      0.02    

+0.02      0.05    

+   0.02      0.08    

+    

+  NEWAXES  CARTES  I  AN  

+0.02   0.0  0.0  0.0  0.0  0.0 

+  NODES  

+0.02   0.02  

+  0.02   0.08  , 

+  0.04   0.02  , 

+  0.04   0.05 

+  0 .  04  0 .  08  , 

+ 

+ 
+ 


Although  this  data  file  Is  correct,  It  Is  quite  evident  that  it 
has  no  specific  usefulness  with  respect  to  this  example. 
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SCALE 
EX  EY  EZ 

If  necessary,  this  entry  defines  scale  factors  (EX,  EY,  EZ)  which 
are  applied  to  the  node  coordinates  Included  in  the  following 
NODES  entry(ies),  during  the  description  of  the  mesh.  To  switch 
off  this  scaling,  the  SCALE  entry  has  to  be  used  with  the 
parameters  0,  0,  0  or  1,  1,  1. 


Test  example.     + 

+   NODES 

+0.00      0.02    

+0.00      0.05    

+   0.00     0.08 

+  0.01      0.02    , 

+0.01      0.08    

+0.02     0.02    

+0.02     0.05    

+  0.02     0.08    , 

+    

+  NEWAXES  CARTES  I  AN  , 

+0.02      0.0   0.0   0.0   0.0   0.0 

+   SCALE 

+   2.0   1.0   1.0    

+  NODES    

+  0.01      0.02    

+0.01      0.08    

+  0.02     0.02    

+0.02     0.05    

+0.02     0.08    

+    

+ 
+ 


This  example  is  a  simple  combination  of  the  entries  NEWAXES  and 
SCALE.  It  allows  the  modification  of  the  thickness  of  the 
electrode  by  simply  modifying  the  scale  factor  (2.00  in  the  above 
file). 


1 1 
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6   ELEMENTS 

RAAAAA  MATERA  NGA 

N11  N12  N13  .    .        .  N1P 

N21  N22  N23  .    .    .    .  N2P 

NK1  NK2  NK3  .    .        .  NKP 

DUMMY  

RBBBBB  MATERB  NGB 

N11  N12  N13 N1Q 

N21  N22  N23 N2Q 

NL1  NL2  NL3 NLQ 

DUMMY  


DUMMY  . 

RZZZZZ  MATERZ  NGZ 

N11  N12  N13  . 

N1R 

N21  N22  N23  . 

N2R 

NM1  NM2  NM3  . 

NMR 

DUMMY 

DUMMY 

This  entry  provides  information  about  the  elements.  It  is 
compulsory.  Its  parameters  are  defined  set  by  set,  each  set 
corresponding  to  successive  elements  of  the  same  type.  The 
various  types  of  elements,  their  topologies  and  their 
corresponding  material  and  geometrical  properties  are  described  in 
chapter  4. 

For  each  set  of  parameters,  the  first  line  contains  the  name  of 
the  element  type,  corresponding  to  RAAAAA,  RBBBBB,  ...,  RZZZZZ  in 
the  above  entry  example,  the  name  of  the  corresponding  material 
data  set,  denoted  MATERx  (see  the  entry  MATERIALS  for  more 
details),  and  the  number  of  the  corresponding  geometrical  property 
set,  denoted  NGx  (see  the  entry  GEOMETRY  for  more  details).  Then, 
each  following  line  contains  the  actual  numbers  of  the  element 
nodes,  given  in  a  compulsory  order  which  corresponds  ta  the 
element  topology  (one  or  two  lines  per  element,  depending  upon  the 
number  of  nodes  belonging  to  it,  the  first  line  being  terminated 
by  the  character  &  if  a  second  line  is  needed).  Finally,  the  set 
has  to  be  closed  by  a  dummy  line.  Moreover,  the  whole  entry 
parameter  list  has  also  to  be  closed  by  a  dummy  line.  Thus,  the 
two  last  lines  of  this  entry  are  always  two  dummy  lines. 
From  the  computational  point  of  view,  elements  are  assembled  in 
the  sequential  order  of  this  entry  parameter  list.  Thus,  in  most 
of  the  cases,  a  great  deal  of  attention  has  to  be  given  to  this 
order,  which  has  an  important  effect  on  the  array  sizes  and  the 
CPU  time  requested  by  the  Job. 
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Test  example.     + 

+  ELEMENTS 


+  R37521  X51  1  

+  1  3682457  

+  

+  R36210  AU4G  2  

+  6811  13  7  9  10  12 

+  

+  

+ 
+ 


The  above  number  I ngs  of  the  nodes  can  be  easily  obtained  using  the 
corresponding  element  descriptions  In  chapter  4.  "X51"  and  "AU4G" 
are  parameters  which  refer  to  the  MATERIALS  entry.  "1"  and  "2"  are 
parameters  which  refer  to  the  GEOMETRY  entry.  Their  meaning  Is 
discussed  in  the  two  following  sections. 
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MATERIALS 

MATER 1 

X11  X12  X13 

MATER2 

X21  X22  X23 


.  X2P 


MATERN 

XN11    XN12   XN13  XN14  XN15   XN16  & 

......  & 

XNM1  XNM2  XNM3  XNM4  XNM5  XNM6  & 
DUMMY  LINE  


This  entry  provides  information  about  the  material  properties.  It 
is  compulsory  and  has  to  be  closed  by  a  dummy  line.  Its  parameters 
are  defined  by  set.  Each  parameter  set  corresponds  to  a  given 
material  and  is  identified  by  a  name,  selected  by  the  user  and 
containing  less  than  8  characters  (MATER1 ,  MATER2,  . . . , MATERN  In 
the  above  example).  This  name  links  a  set  of  elements  from  the 
previous  ELEMENTS  entry  to  the  given  material.  The  list  of  the 
requested  properties  is  provided  in  the  element  description,  In 
chapter  4.  Thus,  as  examples,  this  list  contains  Young's  modulus, 
Poisson's  ratio  and  density  in  the  case  of  a  three-dimensional 
elastic  element,  and  it  contains  the  whole  elastic,  piezoelectric 
and  dielectric  tensors  in  the  case  of  a  piezoelectric  element.  If 
a  line  is  insufficient  to  describe  the  properties  of  a  given 
material,  it  can  be  extended  using  the  character  &.  The  order  in 
which  the  physical  properties  have  to  be  written  Is  compulsory  and 
is  given  below. 

.  elastic  material : 


E  NU  RO 

where  E  is  Young's  modulus, 

densi  ty . 


NU  is  Poisson's  ratio  and  RO  is  the 


piezoelectric  material 


0.0 

0.0 

RO 

0.0 

0.0 

0.0 

S11 

S12 

S13 

S14 

S15 

S16 

S21 

S22 

S23 

S24 

S25 

S26 

S31 

S32 

S33 

S34 

S35 

S36 

S41 

S42 

S43 

S44 

S45 

S46 

S51 

S52 

S53 

S54 

S55 

S56 

S61 

S62 

S63 

S64 

S65 

S66 

d11 

d12 

d13 

d14 

d15 

d16 

d21 

d22 

d23 

d24 

d25 

d26 

d31 

d32 

d33 

d34 

d35 

d36 

EP11 

EP12 

EP13 

0.0 

0.0 

0.0 

EP21 

EP22 

EP23 

0.0 

0.0 

0.0 

EP31 

EP32 

EP33 

0.0 

0.0 

0.0 
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where  [S]  is  the  constant  electric  field  elastic  tensor,  [d]  is 
the  piezoelectric  tensor  and  [EP]  is  the  constant  strain 
dielectric  tensor. 

fluid  material : 

COMP  .0  R0 

where  COMP  is  the  modulus  of  compression  and  RO  is  the  density. 


The  material  property  values  of  classical  or  more  frequently  used 
materials  can  be  written  in  a  special  file,  named  MATER. STD.  If  it 
is  the  case,  these  materials  need  not  be  described  in  the 
MATERIALS  entry  and  their  names  can  be  used  directly  as  parameters 
In  the  ELEMENTS  entry.  If  one  of  the  names  belonging  to  the 
MATER. STD  file  is  used  In  the  MATERIALS  entry,  the  new  property 
values  which  are  listed  as  parameters  will  automatically  replace 
the  old  values  in  the  file  and  will  be  used  by  the  code  during  the 
run.  In  this  case,  a  warning  will  be  printed  in  the  listing  file. 
Finally,  If  a  new  material  name  is  defined  in  the  MATERIALS  entry, 
this  name  and  the  corresponding  property  values  wi  IT  be 
automatically  added  to  the  MATER. STD  file. 

The  MATER. STD  file  can  be  edited  by  the  user  for  any  modification, 
such  as  adding  a  new  material  or  correcting  old  property  values. 
Its  first  line  contains  a  control  integer  (generally  80,  which  is 
the  number  of  columns),  the  date  of  the  last  edition  and  its  total 
number  of  lines.  The  following  lines  contain  the  material  property 
values,  listed  as  defined  above  (1  line  for  elastic  or  fluid 
materials,  13  lines  for  piezoelectric  materials).  Each  of  these 
following  lines  starts  with  two  special  fields,  containing 
respectively  the  material  name  and  the  number  of  lines  set  apart 
for  this  material  (as  examples:  STEEL  1  or  CERAMIC  13  ...).  If 
this  file  is  edited,  the  user  has  to  take  care  of  the  final  total 
number  of  lines  and,  if  needed,  to  correct  its  value  on  the  first 
line  of  the  file.  An  opening  error  can  take  place  during  the 
following  run  of  ATILA  if  this  number  is  wrong. 
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Test  example. 


+ 

+  MATER IALS  

+  AU4G  

+  .714E+11  .344  2780 

+  X51  

+  0.0  0.0  7350 .  0.0  0.0  0.0  

+  0.114E-10  -0.339E-11  -0.410E-11  0.0  0.0  0.0  & 
+  -0.339E-11  0.114E-10  -0.410E-11  0.0  0.0  0.0  & 
+  -0.410E-11  -0.410E-11   0.126E-10  0.0  0.0  0.0  & 

+  0.0  0.0  0.0  0.510E-10  0.0  0.0  &  

+  0.0  0.0  0.0  0.0  0.510E-10  0.0  &  

+  0.0  0.0  0.0  0.0  0.0  0.296E-10  & 

+  0.0  0.0  0.0  0.0-0. 491 E-09  0.0  &  

+  0.0  0.0  0.0  0.491 E-09  0.0  0.0  &  

+  -0.955E-10  -0.955E-10  0.208E-09  0.0  0.0  0.0  & 

+   0.667E-08  0.0  0.0  0.0  0.0  0.0  & 

+   0.0  0.667E-08  0.0  0.0  0.0  0.0  &  

+  0.0  0.0  0.687E-08  0.0  0.0  0.0  , 

+  

+ 
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GEOMETRY 

NG1 

X11  X12 

NG2 

X21  X22  .    .  X2P 


NGM 

XM1  XM2  .       .  XMQ 

DUMMY  LINE  


This  entry  provides  information  about  the  geometrical  properties. 
It  must  be  closed  by  a  dummy  line.  Its  parameters  are  defined  set 
by  set,  each  parameter  set  corresponding  to  a  given  geometry  and 
being  identified  by  a  number  selected  by  the  user  (NG1 , 
NG2,...,NGM  in  the  above  example).  For  a  given  geometry,  the  list 
of  the  requested  properties  is  provided  In  the  element 
description,  in  chapter  4.  Thus,  as  examples,  this  list  contains 
the  thickness  in  the  case  of  a  plate  element,  and  it  contains  the 
radius  of  curvature  in  the  case  of  a  damping  element.  The  order  In 
which  the  geometrical  properties  have  to  be  written  is  compulsory. 


Test  example. 


In  this  case,  the  correct  entry  for  geometrical  properties  is  the 
GEOMETRY  POLARIZATION  entry  which  is  described  In  the  next 
section.  In  fact,  its  format  is  identical  to  the  GEOMETRY  entry 
format,  and  thus  the  example  can  be  considered  as  common  to  the 
two  entr  ies. 
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1- 


GEOMETRY  POLARIZATION  STRING 

NG1 

X11  X12 

NG2 

X21  X22  .    .  X2P 


NGM 

XM1  XM2  .        .  XMQ 

DUMMY  LINE  


This  entry  provides  information  about  the  geometrical  properties 
when  piezoelectric  elements  are  present  in  the  mesh,  and  it  is 
compulsory  In  this  case.  The  GEOMETRY  entry  features  are  all 
relevant  to  this  new  entry  and  must  be  used  for  non-piezoelectric 
elements.  However,  for  piezoelectric  elements,  specific 
parameters  have  to  be  added  to  define  the  polarization  type 
(cartesian,  cylindrical  or  spherical)  and  or lentat Ion. The 
polarization  type  is  specified  by  replacing  STRING  by  CARTESIAN, 
CYLINDRICAL  or  SPHERICAL  in  the  entry  name. 


.  If  STRING-CARTESIAN,  the  polarization  is  homogeneous  in  the 
whole  element  and  the  natural  axes  0X-| ,  0X2  and  OX3  have  to  be 
specified  with  respect  to  the  global  axes  OX,  OY  and  OZ  (OX3  is 
the  polarization  axis  in  the  classical  case  of  a  ceramic 
material).  To  do  this,  the  user  has  to  provide  three  Euler  angles, 
named  ALPHA,  BETA  and  GAMMA,  which  transform  OXYZ  Into  OX3X-,X2  and 
are  defined  as  follows  (see  figure  1): 

-  ALPHA  is  the  angle  of  the  first  rotation,  around  OZ,  which 
transforms  the  OXYZ  system  into  a  OX'Y'Z  system  such  that  OX' 
belongs  to  the  OZX3  plane, 

-  BETA  is  the  angle  of  the  second  rotation,  around  OY',  which 
transforms  the  OX'Y'Z  system  into  a  0X3Y'Z', 

-  GAMMA  is  the  angle  of  the  third  rotation,  around  OX3,  which 
transforms  the  0X3Y'Z'  system  into  the  OX1X2X3  system.  Its  value 
is  dummy  in  the  particular  case  of  a  ceramic  material  (6mm 
symmetry  class) . 

The  entry  parameters  are  ALPHA,  BETA,  GAMMA,  listed  In  this 
compulsory  order.  Attention  has  to  be  given  to  the  exact  meaning 
of  this  "cartesian"  polarization.  In  fact,  If  the  polarization  is 
homogeneous  in  a  three-dimensional  element,  it  is  physically 
homogeneous  in  the  corresponding  part  of  the  material.  On  the 
other  hand,  if  the  polarization  is  homogeneous  and  its  direction 
orthogonal  to  the  symmetry  axis  in  the  case  of  an  ax  isymmetr  ical 
element,  this  polarization  is  physically  cylindrical  In  the 
corresponding  part  of  the  material. 
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.If  STRING-CYLINDRICAL,  the  polarization  is  radial  around  a  given 
axis  0*Z*  in  the  whole  element.  The  user  has  first  to  specify  a 
O'X'Y'Z'  system  which  contains  the  O'Z'  axis.  To  do  this,  he 
provides  the  0'  coordinates  (XO,  YO,  ZO)  as  well  as  the  Euler 
angles  (ALPHA,  BETA,  GAMMA)  which  transform  the  OXYZ  system  into 
the  O'X'Y'Z'  system  and  are  defined  as  follows  (see  figure  2a): 

-  ALPHA  is  the  angle  of  the  first  rotation,  around  OZ,  which 
transforms  the  OXYZ  system  Into  a  OxyZ  system  such  that  Ox 
belongs  to  the  OZX'  plane, 

-  BETA  is  the  angle  of  the  second  rotation,  around  Oy,  which 
transforms  the  OxyZ  system  into  a  OX'yz, 

-  GAMMA  is  the  angle  of  the  third  rotation,  around  OX',  which 
transforms  the  OX'yz  system  into  the  OX'Y'Z'  system. 

It  is  self  evident  that  the  choice  of  the  O'X'  (or  O'Y')  axis  Is 
arbitrary  and  has  to  be  the  simplest.  With  these  new  reference 
axes,  the  angle  of  the  rotation  around  O'Z'  which  brings  the  O'X' 
axis  to  coincidence  with  the  natural  MX3  axis,  for  each  point  M 
which  is  concerned,  is  directly  computed  by  the  code  (see  figure 
2b).  Finally,  the  user  must  also  provide  the  angle  DELTA  which 
brings,  for  every  point  M,  the  MZ'  axis  into  coincidence  with  the 
MX2  natural  axis  (see  figure  2c).  The  value  of  DELTA  is  dummy  in 
the  case  of  a  ceramic  material  (6mm  symmetry  class).  The  entry 
parameters  are  ALPHA,  BETA,  GAMMA,  XO,  YO,  ZO,  DELTA,  listed  In 
this  compulsory  order.  Attention  has  to  be  given  to  the  exact 
meaning  of  this  "cylindrical"  polarization.  In  fact,  if  the 
polarization  is  cylindrical  in  a  three-dimensional  element,  it  is 
physically  cylindrical  in  the  corresponding  part  of  the  material. 
On  the  other  hand,  if  the  polarization  is  cylindrical  around  an 
axis  normal  to  the  plane  of  an  ax isymmetr leal  element  and 
intersecting  the  symmetry  axis  at  the  point  0,  this  polarization 
Is  physically  spherical  around  0  in  the  corresponding  part  of  the 
mater  ial . 


.  If  STRING-SPHERICAL,  the  polarization  is  radial  around  a  given 
point  0'  in  the  whole  element.  The  user  has  first  to  specify  a 
O'X'Y'Z'  system  centered  on  the  origin  0'.  To  do  this,  he  provides 
the  0'  coordinates  (XO,  YO,  ZO)  and  the  Euler  angles  which 
transform  the  OXYZ  system  into  the  O'X'Y'Z'  system  (ALPHA,  BETA, 
GAMMA)  and  are  defined  as  for  the  cylindrical  case  (see  figure 
3a).  It  is  self  evident  that  the  choice  of  these  new  axes  is 
arbitrary  and  has  to  be  the  simplest.  With  these  new  reference 
axes,  the  angles  of  the  rotations  which  bring  the  O'X'  axis  Into 
coincidence  with  the  natural  MX3  axis,  for  each  point  M  which  is 
concerned,  are  computed  by  the  code  (see  figure  2b).  Finally,  the 
user  has  also  to  provide  the  angle  DELTA  which  brings,  for  every 
point  M,  the  MZ'  axis  into  coincidence  with  the  natural  MX2  axis 
(see  figure  2c).  The  value  of  DELTA  is  dummy  in  the  case  of  a 
ceramic  material  (6mm  symmetry  class).  The  entry  parameters  are 
ALPHA,  BETA,  GAMMA,  XO,  YO,  ZO,  DELTA,  listed  in  this  compulsory 
order . 
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Test  example. 


+ 

+  GEOMETRY  POLARIZATION  CARTESIAN 

+  1  , 

+  0.0  0.0  0.0  0.0  0.0  0.0  0.0  

+  2  

+  1 

+  

+ 
+ 


As  can  be  seen  after  "1",  the  three  Euler  angles  are  equal  to- 
zero,  since  the  polarization  direction  and  the  global  OX  axis  are 
identical.  Moreover,  the  value  of  the  thickness  after  "2"  Is 
dummy,  since  the  model  is  ax isymmetr leal . 
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FIGURE  1.  "CARTESIAN"  POLARIZATION. 


a/  definition  of  the  Euler  angles. 


I  * 


*v  ■>  «=> 


b/  examples. 


Ax  I symmetr ica I  case 


Three-dimensional  case. 
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FIGURE  2.  "CYLINDRICAL"  POLARIZATION. 


a/  definition  of  the  geometry. 


^v/ 


b/  definition  of  the  rotation  processed  by  the  code. 


-  rl 


r 


c/  definition  of  the  DELTA  angle. 


22 


Chapter  3.  Description  of  a  data  file 


FIGURE  3.  "SPHERICAL"  POLARIZATION. 


a/  definition  of  the  geometry, 


I*' 


b/  definition  of  the  rotations  processed  by  the  code. 


c/  definition  of  the  DELTA  angle. 


o' 
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REDUCTIONS 
I  N  MDF 


This  entry  defines  the  frontal 
freedom  are  requested  (modal  and 
it  is  compulsory. 


reduction  when  master  degrees  of 
harmonic  analysis).  In  this  case, 


.  I  is  a  control 
provide  NLO  ful I 
backsubst  i tut  ion 
eigenf requencies 
integer  value  of 
integer  N  is 
modes,  which  are 
backsubst  i tut  ion. 
backsubst  i  tut  ion 
paid  to  the  fact 
backsubst I tuted, 


integer.  If  1-0,  backsubst I  tut  Ion 
eigenvectors  from  the  e 
is  processed  to  provide 
of  which  are  above  N 
a  threshold  frequency, 


is  processed  to 

genmode  number  N.  If  1-1, 

NLO  full  eigenvectors  the 

(in  this  case,  N  is  an 

given  in  Hz).  The  control 

mainly  devoted  to  the  elimination  of  the  rigid  body 

provided  by  the  ei genmode  computation,  from  the 

In  some  cases,  it  can  be  used  to  restrict  the 

to  a  specific  given  mode.  Attention  has  to  be 

that  only  the  twenty  first  eigenvectors  can  be 

including  the  rigid  body  modes. 


.  MDF  is  the  chosen  number  of  master  degrees  of  freedom  which  are 
automatically  selected  by  the  code.  Attention  has  to  be  paid  to 
the  fact  that  electrical  potential  degrees  of  freedom  are  always 
condensed,  except  the  applied  potential  in  the  case  of  an  harmonic 
analysis  which  is  a  master  degree  forced  by  the  code.  Moreover, 
manually  selected  degrees  of  freedom  can  be  defined  in  the  set  of 
data  which  specifies  the  boundary  conditions.  Finally,  the  degrees 
of  freedom  corresponding  to  damping  elements  are  all  master 
degrees  and  are  forced  by  the  code.  Then,  the  sum  of  the  total 
master  degree  number  IDF  (automatic  masters  +  applied  potential 
degree  +  degrees  from  damping  elements)  and  the  automatically 
reduced  electrical  degree  number  must  be  lower  than  the  total 
degree  of  freedom  number  IDT. 


Test  example. 


+  REDUCTIONS 
+01  23  

+ 
+ 


Thus,  all  the  available  degrees  of  freedom  are  master  degrees  and, 
due  to  the  fact  that  no  rigid  body  mode  exists  in  this  case.,  the 
five  first  full  eigenvectors  will  be  provided. 
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FREQUENCIES 

F1  F2  F3  FN 

This  entry  is  only  requested  for  an  harmonic  analysis  but  then  is 
compulsory.  Its  parameters  are  the  values  of  the  frequency  for 
which  a  computation  is  needed.  The  maximum  number  of  these  values 
is  20  if  there  is  no  damping  (in-air  harmonic  analysis)  and  10  if 
damping  is  taken  Into  account  (in-water  harmonic  analysis).  The 
parameter  line  can  be  extended  using  the  character  &. 


Test  example. 

There  is  no  test  example  for  this  entry  since  this  entry  is  only 
concerned  with  an  harmonic  analysis. 


RADIATION  MONOPOLAR  or  DIPOLAR 

This  entry  is  only  requested  for  the  modelling  of  a  radiating 
structure  but  then  is  compulsory. 


Test  example. 

There  is  no  test  example  for  this  entry  since  this  entry  Is  only 
concerned  with  an  harmonic  analysis  of  a  radiating  structure. 


PR  I  NT  I NG 
NPR 

This  entry  defines  the  printing  level  NPR,  which  must  equal  to  1, 
2  or  3,  and  corresponds  to  the  amount  of  information  available  on 
the  result  listing.  If  NPR-1 ,  only  the  computation  results  are 
printed.  If  NPR-2,  the  more  important  data  are  recalled  and 
printed  together  with  the  computation  results.  This  level  is  the 
normal  (default)  level.  NPR-3  corresponds  to  the  trace  level  and 
is  only  used  to  detect  faults. 


Test  example. 

For  this  example,  the  default  case  Is  the  best  choice. 
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e   END 

This  entry  marks  the  end  of  the  entry  list  and  is  compulsory, 


Test  example. 


+ 

+  END 

+ 
+ 
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3.4  THE  LOADING  AND  BOUNDARY  CONDITION  DATA 


The  data  which  are  associated  with  the  loading  and  boundary 
conditions  and  which  are  described  hereafter  have  to  be  written  in 
fixed  formats.  Work  is  In  progress  to  provide  these  conditions 
using  specific  entries  in  the  next  ATILA  version.  In  all  the  data 
lines,  a  given  degree  of  freedom  is  marked  by  the  number  D,  with 
the  following  correspondence:  1  for  UX  (or  P),  2  for  UY,  3  for  UZ , 
4  for  NX  (or  V),  5  for  NY  and  6  for  NZ.  Moreover,  lines  which  are 
parallel  to  a  global  coordinate  axis  and  planes  which  are  normal 
to  a  global  coordinate  axis  are  also  marked  by  a  number ,  denoted 
P.  Then,  the  correspondence  is:  1  for  a  plane  normal  to  the  OX 
axis,  2  for  a  plane  normal  to  the  OY  axis,  3  for  a  plane  normal  to 
the  OZ  axis,  4  for  a  line  parallel  to  the  OX  axis,  5  for  a  line 
parallel  to  the  OY  axis  and  6  for  a  line  parallel  to  the  OZ  axis. 


3.4.1  DESCRIPTION  OF  THE  LOADING  CONDITION  DATA. 


In  the  case  of  a  static  analysis,  the  lumped  applied  forces  (or. 
moments)  have  to  be  described  using  one  of  the  two  following 
procedures.  The  procedure  is  selected  using  a  control  number, 
named  CONT,  which  has  to  be  written  first,  In  F10.3  format: 


.CONT, 


Then,  if  CONT-7777.0,  the  nodes  where  the  forces  (or  moments)  are 
applied  are  defined  by  their  coordinates  and,  for  each  applied 
force  (or  moment),  the  following  line  has  to  be  written  (3F10.3, 
2I5,  F10.3): 

X Y Z D..NLO F 


X,  Y,  Z  are  the  coordinates  of  the  fulcrum,  D  is  the  application 
direction  (1  to  3  for  an  applied  force,  4  to  6  for  an  applied 
moment,  following  the  above  mentioned  correspondence  rules),  NLO 
is  the  loading  case  number  (from  1  to  the  value  defined  previously 
bythe  NLOAD  entry),  F  is  the  force  (or  moment)  value. 

If  CONT-8888.0,  the  nodes  where  the  forces  (or  moments)  are 
applied  are  defined  by  their  number  and,  for  each  applied  force 
(or  moment),  the  following  line  has  to  be  issued  (3I5,  F10.3): 

.NODE D..NLO F 


NODE  Is  the  application  node  number,  D,  NLO  and  F  have  the  same 
meanings  as  in  the  previous  case. 

In  each  case,  the  set  of  loading  condition  data  lines  has  to  be 
closed  by  a  dummy  line. 
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WARNING 


If  the  model  is  axlsymmetr leal ,  the  force 
(or  moment)  actual  value  has  to  be  divided 
by  2*P I . 

In  the  case  of  a  distributed  loading,  the 
lumped  nodal  forces  (or  moments)  have  to  be 
computed  using  the  corresponding  element 
shape  functions. 


Test  example. 

The  test  example 
loading  condition 


s  a  modal  analysis  and  does  not  require  any 


3.4.2  DESCRIPTION  OF  THE  BOUNDARY  CONDITION  DATA. 


The  boundary  conditions  are  generally  associated  with 
mechanical  clampings,  symmetry  planes,  electrodes,  pressure 
released  surfaces...  They  are  defined  by  specific  data  lines,  the 
classical  format  of  which  is  (315,  4X,  A1 ) .  The  available  cases 
are  arranged  in  two  groups.  The  first  group  contains  only  boundary 
conditions  which  refer  to  the  global  coordinate  axes.  The  second 
group  contains  boundary  conditions  which  require  specific  local 
coordinates. 


Boundary  conditions  defined  In  global  axes. 

*  The  degree  of  freedom  D  is  deleted  at  the  node  N. 
N D 

*  The  degree  of  freedom  D  is  deleted  for  all  the  nodes  which 
belong  to  a  specific  plane  or  line  P  containing  the  node  N. 

...-N D P 

Thjs  condition  can  be  restricted  to  the  solid  or  the  fluid  part  of 
a  mesh  In  the  case  of  coupled  problems  by  specifying  S  or  F  in  the 
fourth  field  and  using  one  of  the  two  following  lines: 


...-N. 
...-N, 


*  The  degree  of  freedom  D  is  deleted  for  all  the  nodes  in  the 
mesh. 
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*  Degrees  of  freedom  of  type  D  are  Identical  for  all  the  nodes 
which  belong  to  a  specific  line  or  plane  P  containing  the  node  N. 

*  The  degree  of  freedom  D  at  node  N  is  a  master  degree  of  freedom. 


•  *  *  •  PI  •  ■  •  ~U  < 


*  Degrees  of  freedom  of  type  D  are  master  degrees  of  freedom  for 
all  the  nodes  which  belong  to  a  specific  line  or  plane  P 
containing  the  node  N. 


-N...-D 


*  Degrees  of  freedom  of  type  D  are  master  degrees  of  freedom  and 
identical  for  all  the  nodes  which  belong  to  a  specific  line  or 
plane  P  containing  the  node  N. 


.-N. . .-D. . .-P 


*  In  an  harmonic  analysis  of  a  piezoelectric  structure,  the  degree 
of  freedom  associated  with  the  applied  potential  has  to  be  a 
master  degree  of  freedom  and  specifically  marked,  since-  its 
numerical  processing  requires  a  special  storage.  This  selection  is 
realised  using  the  following  data  line: 


where  N  is  a  node  which  belongs  to  the  related  electrode. 
Moreover,  this  marked  degree  of  freedom  is  the  same  for  the  whole 
electrode.  Thus,  If  this  electrode  Is  a  specific  line  or  plane  P, 
the  data  line  can  be  written: 


•  •  •  "**IM  •  •  •  ^"^  •  •  *  ^"r  i 


*  In  an  harmonic  analysis  of  a  piezoelectric  structure  using  the 
mixed  plane  wave-finite  element  method,  the  degree  of  freedom 
associated  with  the  displacement  D  normal  to  the  Interface  between 
the  two  domains  has  to  be  a  master  degree  of  freedom  and 
specifically  marked,  since  its  numerical  processing  requires  a 
special  storage.  This  selection  Is  realised  using  the  following 
data  I ine: 

.  ~.N...-D ;..! 


where  N  is  the  node  which  corresponds  to  D.  Moreover,  this  marked 
degree  of  freedom  is  the  same  for  the  whole  Interface.  Thus,  if 
this  interface  is  a  specific  line  or  plane  P,  the  data  line  can  be 
wr I tten: 
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*  The  degree  of  freedom  of  type  D  at  node  N  Is 
degree  of  freedom  of  the  same  type  D  at  node  M. 


identical  to  the 


N, 


M 


Boundary  conditions  defined  In  local  axes. 


Local  axes  can  be  requested  in  two  different  cases.  Firstly, 
the  displacement  which  is  constrained  can  be  In  a  direction  which 


is  not  parallel  to  the  g 
boundary  condition  has  to 
parallel  to  a  global  axis  or 
to  a  global  axis  (figure 
generally  merged  (figure  3). 


lobal  axis  (figure  1).  Secondly,  a 

be  applied  on  a  line  which  is  not 

on  a  plane  which  is  not  perpendicular 

2).   In  fact,  these  two  cases  are 


\ 


N, 


Figure  1 . 


Figure  2. 


Figure  3. 
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*  To  define  local  axes  at  a  given  node  N,  the  user  has  to  issue 
the  following  lines  (formats  215  and  6F10.3). 


..N...10 


A11 A12 A13 A21 A22 / 

. . .A23 

A11,  A12,  A13  are  the  direction  cosines  of  the  first  local  axis 
and  A21 ,  A22,  A23  those  of  the  second  local  axis.  Then,  degrees  of 
freedom  at  node  N  are  automatically  defined  by  the  code  In  the  new 
local  axes  and  they  can  be  constrained  using  the  same  data  lines 
as  for  the  first  group.  If  the  same  global  axes  have  to  be  defined 
for  all  the  nodes  which  belong  to  the  same  line  or  plane  P 
containing  the  node  N,  the  proceeding  data  set  has  to  be  modified 
by  simply  substituting  for  the  first  line  the  following  line: 


.-N...10. 


*  To   constrain   nodes   belonging   to  a   plane  which   is  not 

perpendicular  to  a  global  axis,  the  corresponding  line  of  the 

first  group  has  to  be  modified  by  substituting  9001  to  P  and 

adding   Immediately  after  a  second   line  which  contains"  the 
direction  cosines  A1 ,  A2  and  A3  of  the  normal  to  this  new  plane 

(format  3F10.3).  Thus,  If  the  degrees  of  freedom  of  type  D  of  all 

the  nodes  which  belong  to  this  plane  have  to  be  deleted,  the  two 
data  I  ines  are: 

. .  .-N D.9001 

A1 A2 A3 


To  constrain  nodes  belonging  to  a  line  which  is  not  parallel  to  a 
global  axis,  the  corresponding  line  of  the  first  group  has  to  be 
modified  by  substituting  9004  to  P  and  adding  Immediately  after  a 
second  line  which  contains  the  direction  cosines  A1 ,  A2  and  A3  of 
this  new  line  (format  3F10.3).  The  proceeding  example  can  be 
immediatly  transposed. 

*  Finally,  the  two  preceeding  cases  can  be  merged  simply  by  mixing 
the  corresponding  data  lines.  Thus,  If  the  preceeding  local  axes 
have  to  be  defined  for  all  the  nodes  belonging  to  the  preceeding 
plane,  the  following  data  lines  have  to  be  Issued: 


.-N... 10. 9001 


.A11 A12 A13 A21 A22 / 

. . .A23 

.  .A1 A2 A3 
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General  comments. 

.  This  data  line  set  has  to  be  closed  by  a  dummy  line. 
.  Constraints  on  different  degrees  of  freedom  can  be  merged  on  the 
same  data  line  by  placing  side  by  side  the  corresponding  integers 
in  the  definition  of  the  D  integer.  Thus,  If  the  same  constraint 
has  to  be  applied  to  UX,  UY  and  NZ,  D  has  to  be  equal  to  126. 
.  Null  flux  condition  is  the  natural  condition  of  the  finite 
element  formulation  used  in  AT  I  LA  and  is  implicit.  This  Is  the 
case  for  the  fluid  velocity  in  fluid  domains  and  for  the  electric 
field  in  piezoelectric  domains. 


Test  example. 


+ 

+  .. 3 

+  .  ..-1  ..  .14 5 

+  ...-6 4 5 

+ 

+ 
+ 


The  first  line  deletes  the  UZ  component  of  the  displacement  field 
(compulsory,  due  to  the  fact  than  LCPDDC-7) .  The  second  line  is 
the  symmetry  condition  along  the  OY  axis.  The  third  line  is  the 
open  electrode  condition  on  the  upper  side  of  the  ceramic  ring. 
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3.5  THE  MIXED  MODEL  SPECIFIC  DATA 


When  a  mixed  finite  element-plane  wave  model  is  used,  the 
finite  element  part  is  described  simply  following  the  guide  lines 
given  in  the  other  sections.  This  section  is  devoted  to  the 
description  of  the  plane  wave  part  of  the  structure,  which  is 
assumed  to  be  organised  as  shown  below. 


///////////. 
///////////. 
///////////. 
/  FINITE  //. 
/  ELEMENT  /. 
/  DOMAIN  //. 
///////////. 
///////////. 
///////////. 


PART  1 


PART  2 


PART  3 


P1 


P2 


P3 


The  plane  wave  part  is  limited  on  one  side  by  its  interface  P1 
with  the  finite  element  domain.  This  interface  has  to  be  a  plane 
and  its  normal  displacement  is  assumed  to  be  uniform.  It  is 
limited  on  the  other  side  by  another  plane  face  P3,  the  normal 
d  i  sp  I  acement  of  wh  i  ch  is  f  ree  or  c  I  amped .  Finally,  it  can  be 
composed  of  one,  two  or  three  parts.  Part  1  is  the  active  part  and 
its  presence  is  compulsory.  Part  3  is  passive  and  in  series  with 
part  1.  It  can  be  missing.  Part  2  is  also  passive  and  in  parallel 
with  part  1.  It  can  be  missing,  but  its  presence  requires  the 
presence  of  part  3.  In  the  classical  case  of  Tonpilz  transducers, 
part  1  is  the  ceramic  stack,  part  2  the  prestress  rod  and  part  3 
the  tailmass,  but  other  schemes  can  be  described  if  they  follow 
the  previous  general  organisation.  The  description  of  this 
structure  is  made  in  the  compulsory  order:  part  1,  part  2,  part  3. 
All  the  constitutive  elements  are  assumed  to  be  of  cylindrical*  or 
conical  shape. 


The  data  lines  are  the  followings.  The  first 
index,  named  IND  (format:  15). 


ine  contains  an 


..IND. 
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If  IND-i,  the  structure  Is  ax Isymmetr leal .  If  IND-2,  the  structure 
is  three-dimensional  with  one  symmetry  plane.  If  IND-3,  the 
structure  is  three-dimensional  with  two  symmetry  planes.  The 
following  lines  describe  the  successive  slices  of  part  1,  2  and  3. 
Their  format  Is  (2I5,  6E10.3).  They  have  to  be  written  under  the 
form: 

. .  ITO.  .  ITH TH 0R1 IR1 0R2 / 

...IR2 PHI 


where: 

ITO  -  NM,  where  N  and  M  are  two  integers 

defined  hereafter, 

ITH  -  0  or  180  and  allows  to  reverse  the 

ceramic  polarization  direction.  Its  value  Is 

dummy  for  other  materials, 

TH  is  the  thickness  of  the  slice, 

0R1  is  the  outside  radius  of  the  slice  first 

face, 

IR1  Is  the  Inside  radius  of  the  slice  first 

face, 

OR2  is  the  outside  radius  of  the  slice  second 

face, 

IR2  is  the  inside  radius  of  the  slice  second 

face, 

PHI  is  the  applied  potential  (beware  of  the 

sign) ,  if  any. 

N  defines  the  type  of  the  slice: 

N  -  1  for  an  elastic  cylinder  (part  1  or  3), 

-  2  for   a   piezoelectric  ceramic  cylinder 
(part  1), 

■  3  for  an  elastic  cone  (part  1  or  3), 

-  -1 ,  -2  or  -3  If  this  slice  Is  the  last 
si  Ice  in  part  3, 

-  4  for  an  elastic  cylinder  (part  2), 

-  5  for  a  conical  cylinder  (part  2), 

-  -4  or  -5  If  this  slice  is  the  last  slice  In 
part  2, 

M  defines  the  slice  material**: 


M  -  1  for  25CD4  steel  or  X5105  ceramic, 

-  2  for  AU4G  or  X9  ceramic, 

-  3  for  brass  (UZ39PB2)  or  X31  ceramic, 

-  4  for   "light"   steel   (RO/100)   or 
ceramic. 


X37 


To  close  this  set,  the  user  has  to  write  down  the  following  data 
I ine  (2I5): 

0. . IP3 
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If  IP3  is  equal  to  zero,  the  normal  displacement  to  P3  Is  free, 
while  if  IP3  is  equal  to  1  this  normal  displacement  is  clamped. 


WARNING:  With  this  version,  only  39  successive  slices  can  be 
accomodated  by  the  code. 


*.  Cylindrical  elements  can  be  replaced  by  elements  with 
different  cross  section  shapes  (squares,  rectangles  ...),  if  the 
plane  wave  assumption  is  still  correct  in  these  cases.  Then,  the 
inside  and  outside  radius  of  the  corresponding  cylindrical  slices, 
as  defined  in  page  34,  have  to  be  computed  by  the  user  to  ensure 
that  the  cross  section  area  values  are  the  actual  values. 

**,.  Material  property  values  are  provided  by  two  .  DATA 
instructions  in  a  subroutine  named  TAC.  The  user  can  modify  them 
by  editing  the  corresponding  FORTRAN  file.  For  elastic  materials, 
values  are  listed  in  the  DATA/SM/  instruction,  with  two  values  for 
each  material  (density  and  Young's  modulus),  successive  materials 
being  put  in  the  M  ascending  order.  For  piezoelectric  materials, 
values  are  listed  In  the  DATA/SP/  instruction,  with  eleven  values 
for  each  material  (density,  constant  electric  field  elastic 
constants  S33,  S13,  S-|-|,  S-j 2 »  S44,  piezoelectric  constants  d33, 
CI3-1  ,  d-|5  and  constant  strain  dielectric  constants  EP33  AND  EP-j  -j ) . 
Successive  materials  are  also  put  in  the  M  ascending  order. 
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3.6  THE  GRAPHIC  DISPLAY  SPECIFIC  DATA 


With  the  ATILA  code,  graphic  displays  can  be  obtained  using 
different  devices,  mainly  plotters  and  workstations.  At  the 
present  time,  results  have  been  obtained  with  BENSON  plotting 
tables,  VERSATEC  Xerox  tables,  Tektronix  4014,  4114,  4115,  as  well 
as  MASSCOMP  and  MICROVAX  computers.  When  plotters  have  to  be  used, 
the  ATILA  run  generates  a  file  which  contains  the  required 
information  for  the  plotter  and  this  plotter  is  generally 
activated  in  a  following  Job.  In  this  case,  graphic  display 
specific  data  have  to  be  added  to  the  data  file  and  are  described 
hereafter.  When  a  workstation  Is  used,  the  graphic  part  of  the  Job 
is  interactive  and  no  data  which  are  related  to  the  graphic 
display  have  to  be  put  in  the  data  file.  Nevertheless,  during  the 
run  of  the  interactive  Job,  questions  are  submitted  to  the  user 
which  correspond  strictly  to  the  following  information. 

Graphic  displays  of  the  structure  in  its  rest  position  or  in 
deformed  shapes  can  be  easily  obtained  using  the  following  set  of 
data.  For  each  graphic  display,  two  data  lines  are  requested  with 
the  respective  formats  (4F10.3,  3I5)  and  (A40,  F10.3).  These  lines 
are: 


.ALPHA BETA XLTH BRTH. .  .NP. .NE1 .  .NE2 

TITLE SCALE 


*  ALPHA  and  BETA  are  angles  (degrees)  which  define  the  projection 
of  the  structure  on  the  drawing  plane.  The  structure  being 
initially  observed  by  the  user  in  the  0Y  direction,  out  of  this 
user,  the  OX  axis  Is  rotated  of  the  ALPHA  angle  around  the  0Z 
axis.  Then,  the  initial  OZ  axis  is  rotated  of  the  BETA  angle 
around  the  original  OX  axis.  Finally,  the  structure  is  projected 
in  the  initial  OXZ  plane.  ALPHA  and  BETA  are  positive  If  the  axis 
are  moving  back  of  the  initial  OXZ  plane. 

*  XLTH  and  BRTH  (in  meters)  define  the  space  occupied  by  the 
graphic  display.  If  BRTH  equals  zero  and  XLTH  is  negative,  ^-XLTH 
is  the  scale  of  the  display  with  respect  to  the  actual  size  of  the 
structure. 


*  NP   defines   the 
correspondences: 


type   of   display,   with   the   following 


NP  -  0  .  Exploded  view  of  the  mesh. 

NP  »  2  .  Exploded  view  of  the  mesh,  with  the 

node  number  ing. 

NP  —2  .  Exploded  view  of  the  mesh,  with  the 

element  numbering. 
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NP  -  3  .  Dashed  view  of  the  mesh,  with  the 
degrees  of  freedom  numbering. 
NP  -  7  .  Simple  display  of   the  mesh   (full 
I Ines) . 

NP  -12  .  Display  of  the  mesh  with  full  outside 
lines  and -dashed  inside  lines. 
NP  -31  to  49  .  Graphic  display  of  the  deformed 
shape  of  the  structure,  in  the  loading  case, 
full  eigenvector  or  excitation  frequency  order. 
The  deformed  shape  is  in  full  line  while  the 
rest  position  is  in  dashed  line.  If  damping  Is 
taken  into  account  in  the  computation,  two 
graphic  displays  have  to  be  issued  for  each 
frequency  (real  and  imaginary  part  of  the 
displacement  field). 

For  NP  -  0,  2,  -2,  3,  7  and  12  the  display  is  provided  by  a 
preprocessor.  If  NP  -  31  to  49,  the  display  is  provided  by  a 
postprocessor.  Several  displays  can  be  simultaneously  requested  in 
each  of  these  two  categories. 

*  NE1  and  NE2  permit  getting  the  graphic  display  of  a  part  of  the 
structure  only.  In  this  case,  NE1  is  the  number  of  the  first  drawn 
element  and  NE2  the  number  of  the  last.  If  these  integers  are 
equal  to  zero,  the  whole  stucture  is  displayed.  Care  must  be  given 
to  the  case  of  fluid  domain,  since  in  this  case  the  finite  element 
variable  Is  the  pressure  which  cannot  be  represented  as  a 
displacement.  So,  fluid  elements  have  to  be  excluded  from  graphic 
displays. 

*  A  specific  TITLE  can  be  given  to  each  graphic  display.  This 
title  cannot  contain  more  than  40  letters. 

*  SCALE  is  a  multiplying  factor  for  the  displacement  amplitude.  It 
allows  adjustment  of  this  amplitude  value  for  sake  of  clarity.  In 
the  case  of  an  eigenmode  computation,  for  which  the  sign  of  the 
whole  displacement  field  is  not  absolutely  defined,  it  allows  also 
to  reversal  of  this  sign  for  the  sake  of  comparison  between 
results  of  different  computations  (  In  the  case  of  piece  part 
analysis,  for  example). 
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CHAPTER  4 
ELEMENT  LIBRARY 


THE  FINITE  ELEMENT  CODE  AT  I  LA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  " INSTITUT  SUPERIEUR  D* ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D* ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN.  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE). 
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4.1  INTRODUCTION 


This  chapter  describes  the  finite  elements  which  are  available 
in  the  ATILA  code  .  For  each  element,  It  provides  Its  present 
name,  its  name  in  the  previous  ATILA  version  (in  parenthesis),  its 
definition,  the  list  of  the  degrees  of  freedom  which  are  attached 
to  its  nodes  (translations,  rotations,  electrical  potential, 
pressure),  its  entry  parameters  (topology,  material  properties, 
geometrical  properties). 

Most  of  the  elements  In  the  ATILA  library  are  Isoparametric. 
Thus,  complex  structures  with  curved  sides  or  faces  can  be  easily 
modelled.  Nevertheless,  best  results  are  obtained  when  these 
elements  have  reasonable  shapes.  Wild  distortions  lead  to 
inaccuracy  because  the  assumptions  made  within  the  code  for 
strains  become  unrealistic.  Extreme  distortion  may  even  "cause 
program  failure.  The  user  should  guard  against  excessive 
distortion  by  adhering  to  the  following  guide  lines  for  2D 
elements  or  for  faces  of  3D  elements: 

-  the  angles  between  adjacent  sides  of  quadrilaterals  should  be 
between  45  and  135  degrees, 

-  the  angles  between  adjacent  sides  of  triangles  should  be  between 
30  and  100  degrees, 

-  In  elements  with  curved  sides,  the  radius  of  curvature  of  each 
side  should  exceed  the  length  of  the  longest  side. 


Warning:  all  data  must  be  given  using  MKS  units  or  in  any  coherent 
unit  system  which  is  deduced  from  the  MKS  system. 


*.  Design  and  test  of  several  other  elements  are  in  progress  in 
the  I. S.E.N.  Acoustics  Laboratory.  Further  information  can  be 
provided  upon  request. 
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4.2  AT  I  LA  ELEMENT  DIRECTORY 


NATURE 

GEOM 

TYPE 

N 

NAME 

COMMENT 

P 

ELAST 1 C 

,  HEXA 

3D 

20 

HEXA20E 

4 

PR  IS 

3D 

15 

PRIS15E 

5 

QUAD 

2D&AX 1 

8 

QUAD08E 

6 

QUAD 

PLATE 

8 

PLAT08E 

7 

QUAD 

SHELL 

8 

FACE08E 

FACET 

8 

TRIA 

2D&AX 1 

6 

TRIA06E 

9 

TRIA 

PLATE 

6 

PLAT06E 

10 

LINE 

SHELL 

3 

SHEL03E 

AXI 

11 

AXI 

5 

TRAN05E 

TRANS 

12 

LINE 

SPRING 

2 

SPRI02E 

13 

PIEZO. 

HEXA 

3D 

20 

HEXA20P 

14 

PRIS 

3D 

15 

PRIS15P 

15 

QUAD 

AXI 

8 

AX  I S08P 

16 

QUAD 

2D 

8 

QUAD08P 

PL.  ST 

17 

TRIA 

AXI 

6" 

AX  I S06P 

18 

TRIA 

2D 

6 

TRIA06P 

PL.  ST 

19 

FLUID 

HEXA 

3D 

20 

HEXA20F 

20 

PRIS 

3D 

15 

PRIS15F 

21 

QUAD 

AXI 

8 

AXIS08F 

22 

TRIA 

AXI 

6 

AXIS06F 

23 

SPECIAL 

QUAD 

3D 

16 

QUAD16I 

INTER 

24 

TRIA 

3D 

12 

TRIA12I 

INTER 

25 

LINE 

AXI 

6 

AXIS06I 

INTER 

26 

QUAD 

3D 

8 

QUAD08M 

MONOP 

27 

TRIA 

3D 

6 

TR I A06M 

MONOP 

28 

LINE 

AXI 

3 

AX  I S03M 

MONOP 

29 

LINE 

AXI 

3 

AX  I S03D 

DIPOL 

30 

GEOM :geometr leal  shape 
of  the  element. 

N   : number  of  nodes. 

NAME: name  of  the  element 
as  required  for  the 
entry. 

P   :page  number. 


HEXA  rhexaedra. 

PRIS  :trlangular  prism. 

QUAD  rquadr I  lateral . 

TRIA  :triangle. 

INTER: Interface  element. 

MONOP monopolar  damper. 

DIPOL:di polar  damper. 

PL.ST:plane  strain. 

TRANS: trans  It  ion  element. 
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4.3  ELEMENT  DESCRIPTION 


Descriptions  of  a  variety  of  elements  are  provided  on  the 
fol lowing  pages. 


Chapter  4.  Element  Ij 


HEXA20E  (R37110). 


DESCRIPTION. 

.HEXA20E   is  a  twenty  node   Isoparametric  hexaedra,  used  for 
homogeneous  isotropic  elastic  materials. 

.The  degrees  of   freedom  are,   for   each  node:   ux,  uy,  uz 
(3  translat  ions) . 

TOPOLOGY. 
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19   . 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

HEXA20E  MATER 

N1    N2   N3  N4  N5  N6   N7  N8  N9  N10  N11    N12   N13  N14  N15  N16  N17  N18  N19 

N20 

where  N1....N20  are  the  actual  node  numbers  corresponding  to  the 
integers  1....20  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3) . 

."GEOMETRY"  entry  parameters. 
No. 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3 
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PRIS15E  (R37210). 


DESCRIPTION. 

.PRIS15E  Is  a  fifteen  node  Isoparametric  triangular  prism,  used 
for  homogeneous  isotropic  elastic  material. 


.The  degrees  of   freedom  are,   for  each 
(3  translat  ions) . 


node 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

PRIS15E  MATER 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12  N13  N14  N15 

where  N1....N15  are  the  actual  node  numbers  corresponding  to  the 
integers  1....15  on  the  above  figure. 

."-MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3). 

."GEOMETRY"  entry  parameters. 
No. 

COMMENTS . 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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QUAD08E  (R36210). 


DESCRIPTION. 

.QUAD08E  is  an  eight  node  isoparametric  quadrilateral,  used  for 
homogeneous  isotropic  elastic  material  and  for  plane  stress,  plane 
strain  or  ax isymmetr ical  analysis. 

.The  degrees  of  freedom  are,  for  each,  node:  ux,  uy 
(2  translations). 

TOPOLOGY. 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD08E  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3) . 

."GEOMETRY"  entry  parameter. 

T 

where  T  is  the  thickness  of  the  element  (ignored  In  the  case  of 

plane  strain  and  ax Isymmetr lea  I  analyses). 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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Chapter  4.  Element  library 


PLAT08E  (R38110). 


DESCRIPTION. 

.PLAT08E  is  an  eight  node  flat  quadrilateral  plate  element,  used 
for  homogeneous  isotropic  elastic  material.  It  relies  on  the 
classical  Love-Kirchhof f  hypothesis. 


.The  degrees  of  freedom  are,   for  each 
(1  translation,  2  rotations). 


node 


u 


2» 


n 


x» 


TOPOLOGY. 


.  .  X 


ENTR I ES . 


."ELEMENTS"  entry  parameters. 

PLAT08E  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

T 

where  T  is  the  thickness  of  the  element. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  longitudinal  and  lateral  dimensions  of  the  whole  plate  have 

to  be  at  least  ten  times  the  thickness. 

.This  element  can  be  superimposed  to  the  QUAD08E  element  (plane 

stress  condition)  to  generate  a  flat  shell  element  (facet). 


Chapter  4.  El ement 


FACE08E  (R44210). 


DESCRIPTION. 

.FACE08E  Is  an  eight  node  quadrilateral  shell  element,  used  for 
homogeneous  Isotropic  elastic  material.  It  relies  on  the  classical 
Love-Ki rchhof f  hypothesis  and  has  to  be  flat  (facet). 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy>  u2,  nx,  ny 
(3  translations,  2  rotations). 

TOPOLOGY. 
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0  .  .  .  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

FACE08E  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 
integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3) . 

."GEOMETRY"  entry  parameter. 

T 

where  T  Is  the  thickness  of  the  element  . 

COMMENTS . 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  longitudinal  and  lateral  dimensions  of  the  whole  shell  have 

to  be  at  least  ten  times  the  thickness. 

.This  element  is  equivalent  to  the  superposition  of  the  QUAD08E 

element  (plane  stress  condition)  and  PLAT08E  element. 
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TRIA06E  (R36110). 


DESCRIPTION. 

.TRIA06E  Is  a  six  node  isoparametric  triangle,  used  for 
homogeneous  isotropic  elastic  material  and  for  plane  stress,  plane 
strain  or  ax isymmetr ical  analysis. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy 
(2  translations). 

TOPOLOGY. 


3 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 
TRIA06E  MATER  NG 
N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

T 

where  T  is  the  thickness  of  the  element  (ignored  in  the  case  of 

plane  strain  and  ax isymmetr ical  analyses). 

COMMENTS . 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 


Chapter  4.  Element 


PLAT06E  (R38401). 


DESCRIPTION. 

.PLAT06E  is  a  six  node  flat  triangular  plate  element,  used  for 
homogeneous  isotropic  elastic  material.  It  relies  on  the  classical 
Love-Ki rchhof f  hypothesis. 

.The  degrees  of  freedom  are,  for  each  node:  uz,  nx,  ny 
(1  translation,  2  rotations). 

TOPOLOGY. 


ENTRIES. 

."ELEMENTS"  entry  parameters. 
PLAT06E  MATER  NG 
N1  N2  N3  N4  N5  N6 

where  N1,...N6  are  the  actual  node  numbers  corresponding  to  the 
integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  Is  Poisson's  ratio  and  RO  Is 

the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

T 

where  T  is  the  thickness  of  the  element. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  longitudinal  and  lateral  dimensions  of  the  whole  plate  have 

to  be  at  least  ten  times  the  thickness. 

.This  element  can  be  superimposed  to  the  TRIA06E  element  (plane 

stress  condition)  to  generate  a  flat  shell  element  (facet). 
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Chapter  4.  Element  library 


SHEL03E  (R42123). 


DESCRIPTION. 

.SHEL03E  Is  a  three  node  ax Isymmetr leal  thin  shell  element,  the 
formulation  of  which  takes  account  of  its  double  curvature.  It 
relies  on  the  classical  Love-Kirchhof f  hypothesis. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uv,  n2 
(2  translations,  1  rotation). 

TOPOLOGY. 


0  .  .  .  X 


ENTRIES. 

."ELEMENTS"entry  parameters. 

SHEL03E  MATER  NG 

N1  N2  N3 

where  N1 ,  N2,  N3  are  the  actual  node  numbers  corresponding  to  the 

integers  1,  2,  3  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

T  R 

where  T  Is  the  thickness  of  the  element  and  R  its  radius  of 

curvature  (around  the  axis). 


COMMENTS . 

.The  longitudinal  and  radial  dimensions  of  the  whole  shell  have  to 
be  at  least  ten  times  the  thickness. 


1 1 


Chapter 


4.  Element! 


TRAN05E  (R43000). 


DESCRIPTION. 

.TRAN05E  is  a  five  node  ax isymmetr ica I  element  which  ensures  the 
transition  between  quadrilateral  or  triangular  ax isymmetr leal 
elements  and  the  previous  thin  shell  ax isymmetr ica I  element.  It 
allows  an  accurate  description  of  their  interaction.  It  includes  a 
solid  part  (nodes  1,  2,  and  3)  and  a  shell  part  (nodes  4  and  5). 


.The  degrees  of  freedom  are,  for  nodes  1,  2  and  3 
(2  translations),   and   for   nodes  4  and  5:   u 
translations,  1  rotation). 


x» 


Ux 


u 
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<! 


TOPOLOGY. 


0  .  . 


4  ...  5 


ENTRIES. 

."ELEMENTS"entry  parameters. 

TRAN05E  MATER  NG 

N1  N2  N3  N4  N5 

where  N1 ,   N2   ,   N3,   N4,   N5  are  the  actual   node  numbers 

corresponding  to  the  integers  1,  2,  3,  4,  5  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  NU  RO 

where  E  is  Young's  modulus  (Pa),  NU  is  Poisson's  ratio  and  RO  is 

the  density  (kg/m3). 


."GEOMETRY"  entry  parameter. 

T  B 

where  T  is  the  thickness  of  the  element  and  R 

curvature  (around  the  axis). 


its  radius  of 


COMMENTS. 

.The  longitudinal  and  radial  dimensions  of  the  whole  shell  have  to 
be  at  least  ten  times  the  thickness. 


12 


Chapter  4.  Element  library 


SPRI02E  (R30100). 


DESCRIPTION. 

.SPRI02E  is  a  two  node  spring  element,  which  is  only  able  to 
transmit  to  Its  two  limiting  nodes  a  restoring  force  due  to  a 
change  of  Its  length.  Its  classical  use  Is  for  the  modelling  of 
prestress  bolts  in  transducer  design. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  uz 
(3  translat  ions) . 

TOPOLOGY. 


.  .  X  1 2 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

SPRI02E  MATER  NG 

N1  N2 

where  N1  and  N2  are  the  actual  node  numbers  corresponding  to  the 

integers  1  and  2  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

E  .0  RO 

where  E  is  Young's  modulus  (Pa)  and  RO  Is  the  density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

S 

where  S  is  the  cross  section  area. 

COMMENTS. 

.In  ax isymmetr leal  models,  the  actual  area  S  has  to  be  divided  by 
2*pi.  In  three-dimensional  models,  if  the  rod  belongs  to  N 
symmetry  planes,  its  actual  area  S  has  to  be  divided  by  N+1 . 
.This  element  is  nothing  more  than  a  classical  spring.  Attention 
has  to  be  paid  to  spurious  modes  which  can  be  due  to  the  simple 
spring  hypothesis. 
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Chapter  4.  Element 


HEXA20P  (R37511). 


DESCRIPTION. 

.HEXA20P  is  a  twenty  node  Isoparametric  hexaedra,  designed  for  any 
piezoelectric  material. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  uz,  v 
(3  translations,  1  electrical  potential). 

TOPOLOGY. 


7  ...  .20  ....  8 
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ENTR I ES . 

."ELEMENTS"  entry  parameters. 

HEXA20P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12  N13  N14  N15  N16  N17  N18  N19 

N20 

where  N1....N20  are  the  actual  node  numbers  corresponding  to  the 
integers  1....20  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elast ic.piezoelectr Ic  and  dielectric  tensors  (see  the 

detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  issued  under  the  special  form 

"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 
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Chapter  4.  Element  library, 


PRIS15P  (R37512). 


DESCRIPTION. 

.PRIS15P   is  a  fifteen  node   isoparametric  triangular  prism, 
designed  for  any  piezoelectric  material. 


.The  degrees  of  freedom  are,  for  each  node:  u 
(3  translations,  1  electrical  potential). 

TOPOLOGY. 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

PRIS15P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12  N13  N14  N15 

where  N1,...N15  are  the  actual  node  numbers  corresponding  to  the 
integers  1....15  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elast lc,p iezoe lectr ic  and  dielectric  tensors  (see  the 

detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  issued  under  the  special  form 

"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 
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Chapter  4.  Element 


AXIS08P  (R37521). 


DESCRIPTION. 

.AXIS08P  is  an  eight  node  isoparametric  quadrilateral,  designed 
for  any  piezoelectric  material  and  used  only  for  ax isymmetr leal 
analysis. 

.The  degrees  of  freedom  are,  for  each  node:  ux,  uy,  v 
(2  translations,  1  electrical  potential). 

TOPOLOGY. 


0  .  .  .  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS08P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1,...N8  are  the  actual  node  numbers  corresponding  to  the 

integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 

detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS, 


.Generally,  the  aspect 
.The  entry  "GEOMETRY" 
"GEOMETRY  POLARIZATION 


ratio  has  to  be  less  than  3. 
has  to  be  issued  under  the  special 
(CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 


form 
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Chapter  4.  Element  library 


QUAD08P  (R37501). 


DESCRIPTION. 


.QUAD08P  is  an  eight  node  isoparametric  quadrilateral,  designed 
for  any  piezoelectric  material  and  used  only  for  plane  strain 
analysis. 


.The  degrees  of   freedom  are,   for   each  node: 
(2  translations,  1  electrical  potential). 

TOPOLOGY. 
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Ux 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD08P  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1 ...  .N8  are  the  actual  node  numbers  corresponding  to  the 

integers  1,...8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 

detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS. 


.Generally,  the  aspect 
.The  entry  "GEOMETRY" 
"GEOMETRY  POLARIZATION 


ratio  has  to  be  less  than  3. 
has  to  be  issued  under  the  sped 
(CARTESIAN,  CYLINDRICAL,  SPHERICAL) 


form 
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Chapter  4.  Element 


AXIS06P  (R37522). 


DESCRIPTION. 

.AXIS06P  is  a  six  node  isoparametric  triangle,  designed  for  any 
piezoelectric  material  and  used  only  for  ax isymmetr leal  analysis. 

.The   degrees  of   freedom  are,   for   each   node:   ux,   uy,   v 
(2  translations,  1  electrical  potential). 

TOPOLOGY. 
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0  .  .  .  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 
AXIS06P  MATER  NG 
N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 
integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 

detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  issued  under  the  special  form 

'GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 
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Chapter  4.  Element  library 


TRIA06P  (R 37502). 


DESCRIPTION. 

.TRIA06P  is  a  six  node  isoparametric  triangle,  designed  for  any 
piezoelectric  material  and  used  only  for  plane  strain  analysis. 

.The  degrees  of   freedom  are,   for   each  node:   ux,   uy,   v 
(2  translations,  1  electrical  potential). 

TOPOLOGY. 


6 


0  .  .  .  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

TRIA06P  MATER  NG 

N1  N2  N3  N4  N5  N6 

where  N1,...N6  are  the  actual  node  numbers  corresponding  to  the 

integers  1,...6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

The  whole  elastic,  piezoelectric  and  dielectric  tensors  (see  the 

detailed  description  of  the  entry:  MATERIALS). 

."GEOMETRY"  entry  parameters. 

The  new  origin  coordinates  and  the  Euler  angles  which  define  the 
natural  axes  of  the  material  (see  the  detailed  description  of  the 
entry:  GEOMETRY  POLARIZATION). 


COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.The  entry  "GEOMETRY"  has  to  be  issued  under  the  special  form 

"GEOMETRY  POLARIZATION  (CARTESIAN,  CYLINDRICAL,  SPHERICAL)". 
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Chapter  4.  Element 


HEXA20F  (R 37810) 


DESCRIPTION. 

.HEXA20F  is  a  twenty  node  isoparametric  hexaedra,  used  for  fluid 
homogeneous  media. 

.The  degree  of  freedom  is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

HEXA20F  MATER 

N1    N2  N3  N4  N5  N6  N7  N8  N9  N10  N11    N12   N13  N14  N15  N16  N17  N18  N19 

N20 

where  N1,...N20  are  the  actual  node  numbers  corresponding  to  the 
integers  1,...20  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  Is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameters. 
No. 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3, 
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PRIS15F  (R37820). 


DESCRIPTION. 

.PRIS15F  is  a  fifteen  node  isoparametric  triangular  prism,  used 
for  fluid  homogeneous  media. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 
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ENTR | ES . 

."ELEMENTS"  entry  parameters. 

PRIS15F  MATER 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12  N13  N14  N15 

where  N1....N15  are  the  actual  node  numbers  corresponding  to  the 
integers  1....15  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  (Pa),  and  RO  is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameters. 
No. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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AXIS08F  (R37531). 


DESCRIPTION. 

.AXIS08F  is  an  eight  node  Isoparametric  quadrilateral,  used  for 
fluid  homogeneous  media  and  for  ax Isymmetr leal  analysis. 

.The  degree  of  freedom  is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


.  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS08F  MATER 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1....N8  are  the  actual  node  numbers  corresponding  to  the 

integers  1,...8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  is  the 

density  (kg/m^). 

."GEOMETRY"  entry  parameters. 
No_ 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 
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AXIS06F  (R37532). 


DESCRIPTION. 

.AXIS06F  is  a  six  node  Isoparametric  triangle  used  for  fluid 
homogeneous  media  and  for  ax isymmetr leal  analysis. 

.The  degree  of  freedom  is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


6 


0  .  .  .  X 


ENTR I ES . 

."ELEMENTS"  entry  parameters. 

AXIS06F  MATER 

N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 

integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameters. 
No^. 

COMMENTS. 


.Generally,  the  aspect  ratio  has  to  be  less  than  3. 


23 


Chapter  4.  Elemenl 


QUAD16I  (R37910). 


DESCRIPTION. 

.QUAD16I  is  a  sixteen  node  isoparametric  quadrilateral  element 
which  has  to  be  used  to  ensure  the  matching  between  solid  and 
fluid  meshes  along  their  interface  in  the  case  of  three- 
dimensional  analyses.  It  includes  8  solid  nodes  and  8  fluid  nodes 
which  have  the  same  coordinates  as  the  solid  nodes.  It  has  no 
th  ickness. 


.The  degrees  of  freedom  are,  for  each  solid  node:  u 
(3  translations),  and  for  each  fluid  node:  p  (the  pressure) 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD16I 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12  N13  N14  N15  N16 

where  N1....N16  are  the  actual  node  numbers  corresponding  to  the 
integers  1....16  on  the  above  figure. 


."MATERIALS"  entry  parameters 
No. 

."GEOMETRY"  entry  parameters. 
No. 


COMMENTS. 

.The  nodes  1,  2,  3,  4,  9,  10,  11,  12  are  fluid  nodes.  The  nodes  5, 
6,  7,  8,  13,  14,  15,  16  are  solid  nodes. 

.Attention  must  be  paid  to  the  orientation  of  the  0Z  axis,  from 
the  fluid  domain  to  the  solid  domain. 
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TRIA12I  (R37920). 


DESCRIPTION. 

.TRIM 2 1  is  a  twelve  node  Isoparametric  triangular  element  which 
has  to  be  used  to  ensure  the  matching  between  solid  and  fluid 
meshes  along  their  Interface  In  the  case  of  three-dimensional 
analyses.  It  includes  6  solid  nodes  and  6  fluid  nodes  which  have 
the  same  coordinates  as  the  solid  nodes.  It  has  no  thickness. 

.The  degrees  of  freedom  are,  for  each  solid  node:  ux,  uY,  uz 
(3  translations),  and  for  each  fluid  node:  p  (the  pressure)  . 

TOPOLOGY. 
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FLUID  MESH 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

TRIM  2 1 

N1  N2  N3  N4  N5  N6  N7  N8  N9  N10  N11  N12 

where  N1....N12  are  the  actual  node  numbers  corresponding  to  the 

integers  1 ... .12  on  the  above  figure. 

."MATERIALS"  entry  parameters. 
No. 

."GEOMETRY"  entry  parameters. 
No. 

COMMENTS . 

.The  nodes  1,  2,  3,  7,  8,  9  are  fluid  nodes.  The  nodes  4,  5,  6, 
10,  11,  12  are  solid  nodes. 

.Attention  must  be  paid  to  the  orientation  of  the  OZ  axis,  from 
the  fluid  domain  to  the  solid  domain. 
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AXIS06I  (R37541). 


DESCRIPTION. 

.AXIS06I  is  a  six  node  isoparametric  element  which  has  to  be  used 
to  ensure  the  matching  between  solid  and  fluid  meshes  along  their 
interface  in  the  case  of  ax Isymmetr leal  analyses.  It  includes  3 
solid  nodes  and  3  fluid  nodes  which  have  the  same  coordinates  as 
the  solid  nodes.  It  has  no  thickness. 


.The  degrees  of  freedom  are,   for  each  solid  node:  u 
(2  translations),  and  for  each  fluid  node:  p  (the  pressure) 
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TOPOLOGY , 
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ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS06I 

N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 

integers  1....6  on  the  above  figure. 

."MATERIALS"  entry  parameters. 
No. 

."GEOMETRY"  entry  parameters. 
No. 

COMMENTS. 

.The  nodes  1,  2,  5  are  fluid  nodes.  The  nodes  2,  3,  6  are  solid 

nodes. 

.Attention  must  be  paid  to  the  orientation  of  the  OY  axis,  from 

the  fluid  domain  to  the  solid  domain. 
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Chapter  4.  Element  library 


QUAD08M  (R 37930). 


DESCRIPTION. 

.QUAD08M  is  an  eight  node  isoparametric  quadrilateral  damping 
element  used  to  take  account  of  a  radiation  condition.  It  is 
restricted  to  the  monopolar  (spherical)  damping  condition.  It  has 
to  be  spherical  and  attached  to  the  outside  surface  of  a  three- 
dimensional  fluid  mesh. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

QUAD08M  MATER  NG 

N1  N2  N3  N4  N5  N6  N7  N8 

where  N1,...N8  are  the  actual  node  numbers  corresponding  to  the 

integers  1....8  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  is  the  radius  of  curvature  of  the  element. 

COMMENTS . 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.With  respect  to  the  orientation  of  the  whole  mesh,  Its  normal  has 

to  point  outwards. 
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TRIA06M  (R 37940). 


DESCRIPTION. 

.TRIA06M  is  a  six  node  isoparametric  triangular  damping  element 
used  to  take  account  of  a  radiation  condition.  It  Is  restricted  to 
the  monopolar  (spherical)  damping  condition.  It  has  to  be 
spherical  and  attached  to  the  outside  surface  of  a  three- 
dimensional  fluid  mesh. 

.The  degree  of  freedom  is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

TRIA06M  MATER  NG 

N1  N2  N3  N4  N5  N6 

where  N1....N6  are  the  actual  node  numbers  corresponding  to  the 

Integers  1,...6  on  the  above  figure. 

."MATERIALS  entry  parameters. 

COMP  .0  R0 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  is  the  radius  of  curvature  of  the  element. 

COMMENTS. 

.Generally,  the  aspect  ratio  has  to  be  less  than  3. 

.With  respect  to  the  orientation  of  the  whole  mesh,  Its  normal  has 

to  point  outwards. 
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Chapter  4.  Element  library, 


AXIS03M  (R37551). 


DESCRIPTION. 

.AXIS03M  is  a  three  node  Isoparametric  damping  element  used  to 
take  account  of  a  radiation  condition.  It  is  restricted  to  the 
monopolar  (spherical)  damping  condition.  It  has  to  be  circular  and 
attached  to  the  outside  surface  of  an  ax isymmetr lea  I  fluid  mesh. 

.The  degree  of  freedom  is,  for  each  node:  p  (the  pressure). 

TOPOLOGY . 


0  .  .  .  X 


ENTRIES. 

."ELEMENTS"  entry  parameters. 

AXIS03M  MATER  NG 

N1  N2  N3 

where  N1 ,  N2  and  N3  are  the  actual  node  numbers  corresponding  to 

the  integers  1,  2  and  3  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  RO 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  is  the  radius  of  curvature  of  the  element. 

COMMENTS . 

.To  ensure  the  correct  orientation  of  this  element,  its  local  OX 
axis  must  be  oriented  in  the  positive  direction  of  the  global  OX 
axis.  Thus,  the  global  x  coordinate  of  node  2  must  be  higher  than 
the  global  x  coordinate  of  node  1  on  the  above  figure. 
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AXIS03D  (T37551). 


DESCRIPTION. 

.AXIS03D  is  a  three  node  isoparametric  damping  element  used  to 
take  account  of  a  radiation  condition.  It  includes  the  dipolar 
damping  condition.  It  has  to  be  circular  and  attached  to  the 
outside  surface  of  an  ax isymmetr leal  fluid  mesh. 

.The  degree  of  freedom  Is,  for  each  node:  p  (the  pressure). 

TOPOLOGY. 


0  .  .  .  X 


ENTRIES. 

."ELEMENTS"entry  parameters. 

AXIS03D  MATER  NG 

N1  N2  N3 

where  N1 ,  N2  and  N3  are  the  actual  node  numbers  corresponding  to 

the  integers  1,  2  and  3  on  the  above  figure. 

."MATERIALS"  entry  parameters. 

COMP  .0  R0 

where  COMP  is  the  modulus  of  compression  (Pa)  and  RO  is  the 

density  (kg/m3). 

."GEOMETRY"  entry  parameter. 

R 

where  R  is  the  radius  of  curvature  of  the  element. 

COMMENTS. 


.To  ensure  the  correct  orientation  of  this  element,  its  local  OX 
axis  must  be  oriented  in  the  positive  direction  of  the  global  OX 
axis.  Thus,  the  global  x  coordinate  of  node  2  must  be  higher  than 
the  global  x  coordinate  of  node  1  on  the  above  figure. 
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CHAPTER  5 
DESCRIPTION  OF  AN  AT  I  LA  RUN 


THE  FINITE  ELEMENT  CODE  AT  I  LA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  "INSTITUT  SUPERIEUR  D* ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D* ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN.  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE. 
FRANCE). 


Chapter  5.  Description  of  an  AULA  run 


5.1  INTRODUCTION 


The  general  organisation  of  an  ATILA  Job  Is  described  In 
section  1.3.  This  chapter  provides  the  corresponding  organisation 
of  an  ATILA  run  in  the  special  case  of  a  VAX  or  MICROVAX  computer 
using  the  VMS  operating  system.  The  same  organisation  is  also 
correct  for  other  computer  types  but  with  different  procedures. 
Moreover,  the  part  of  this  chapter  which  Is  devoted  to  the  graphic 
displays  implies  that  the  GKS  library  Is  available  and  can  be 
linked  during  the  corresponding  part  of  the  run,  the  user  being 
then  working  interactively  with  the  help  of  a  graphic  terminal 
(Tektronix  4014,  VT240  ...). 
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5.2  THE  MAIN  RUN 


To  run  an  ATILA  job,  the  user  must  have  written  a  data  file  as 
described  in  chapter  3.  This  data  file  must  be  named  JOBNAME.ATI, 
ATI  being  the  job  type.  Then  the  user  must  run  the  automatic 
program  generator,  named  PGEN,  using  the  entry: 

RUN  PGEN 

During  this  run,  several  questions  are  asked  by  the  code  which 
allow  the  specification  of  the  Job  name  (JOBNAME)  and  the  job  type 
(ATI),  the  verification  of  the  analysis  type  which  has  been 
identified  by  the  code,  the  requirement  of  special  data  storages 
for  a  job  restart  (only  In  the  case  of  an  harmonic  analysis)  or  a 
postprocessing  (graphic  displays  ...).  PGEN  provides  three  files 
named  JOBNAME. FOR,  J0BNAME.COM  and  JOBNAME. LST.  JOBNAME. FOR 
contains  the  ATILA  main  program.  Generally,  it  can  be  used  wilhout 
any  modification.  Nevertheless,  it  has  to  be  edited  In  few  special 
cases.-  the  modal  analysis  of  an  elastic  structure  flooded  in  an 
incompressible  fluid  (MOD  5),  the  radiated  pressure  computation 
for  a  given  displacement  field  (HAR  3)  and  the  use,  by  advanced 
users,  of  new  algorithms  which  are  currently  under  test  (MOD  1, 
MOD  4,  ...).  In  these  cases,  footnotes  are  available  at  the  end  of 
the  corresponding  sections  in  chapter  2  which  provide  all  the 
necessary  information.  J0BNAME.COM  is  used  to  compile  and  link  the 
program.  It  need  not  be  edited.  Finally,  JOBNAME. LST  contains  a 
copy  of  the  data  file,  as  read  during  the  PGEN  run,  as  well  as 
information  about  the  degree  of  freedom  numbering.  It  allows  a 
very  efficient  check  of  the  data  file  and  mainly  of  the  boundary 
conditions.  Generally,  the  user  must  edit  this  list  and  perform 
the  corresponding  check. 

After  the  PGEN  run,  the  user  has  to  compile  the  main  progam  and 
link  it  to  the  requested  libraries.  This  step  is  simply  realised 
using  the  entry: 

6 JOBNAME 

which  provides  the  executable  (JOBNAME.EXE).  Finally,  the  user  has 
to  run  the  program  using  the  entry: 

RUN  JOBNAME 

During  this  run,  the  code  asks  for  the  Job  name  (JOBNAME)  and  the 
job  type  (ATI)  and  information  appears  on  the  user's  screen  which 
describes  the  current  status  (number  of  the  element  the  code  Is 
currently  assembling,  name  of  the  solver  the  code  Is  currently 
using  ...)  and  provides  the  timing.  Moreover,  the  code  writes  down 
several   flies  named  JOBNAME. LST,   JOBNAME. WRK,   JOBNAME. SY2  and 
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J0BNAME.SY4.  JOBNAME.LST  contains  a  partial  copy  of  the  data  file, 
depending  upon  the  printing  level  selected  by  the  user  in  this 
data  file,  and  all  the  results  (eigenf requencies,  displacement  or 
pressure  fields  associated  with  successive  elgenmodes  or 
excitation  frequencies,  electrical  impedances  ...  ).  The 
JOBNAME.WRK  and  J0BNAME.SY2  files  are  working  arrays  the  code 
always  needs.  They  contain  useful  information  at  the  end  of  the 
run  only  if  a  restart  has  to  be  performed.  The  J0BNAME.SY4  file 
contains  several  arrays  the  code  needs  for  postprocessing. 
Particularly,  it  can  provide  the  PLO  array  which  contains  the 
successive  displacement  and/or  pressure  fields  as  well  as  the 
CPDDC  array  which  contains  the  node  coordinates  and  the  degree  of 
freedom  numbering.  This  file  can  be  easily  read  by  advanced  users, 
by  editing  the  JOBNAME.FOR  file,  deleting  all  the  subroutine  calls 
and  replacing  the  WRITE  entries  devoted  to  the  SY4  file  by  the 
corresponding  READ  entries. 
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5.3  THE  RESTART 


To  restart  an  ATILA  job  is  often  necessary  when  the  user  is 
performing  an  harmonic  analysis,  since  more  than  10  or  20 
different  frequencies  are  often  required  in  the  frequency  range  of 
interest.  The  restart  allows  the  computation  of  the  frequency 
response  using  previously  assembled  matrices  and  often  saves  a 
large  amount  of  CPU  time.  Obviously,  it  can  only  be  performed  if  a 
full  run  has  been  realised  before  and  if  the  required  arrays  have 
been  stored.  Then,  the  process  is  exactly  the  same  as  for  a  main 
run,  except  that  the  user  has  to  select  the  restart  during  the  run 
of  PGEN,  when  the  corresponding  question  appears  on  the  screen. 
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5.4  PRE  AND  POSTPROCESSING 


In  this  version  of  the  ATILA  code,  the  available  pre  and 
postprocessors  are  mainly  devoted  to  the  graphic  display  using  the 
GKS  package,  although  many  specific  postprocessors  (farfield 
computation,  transmitting  voltage  response  and  directivity 
computations  and  displays,  stress  computations  and  displays, 
isovalue  displays  ...)  are  currently  under  test  but  are  not 
documented  for  the  present  time.  At  the  preprocessing  level,  the 
code  can  provide  graphic  displays  of  the  mesh  with  different  types 
of  numbering  (nodes,  elements,  degrees  of  freedom).  At  the 
postprocessing  level,  it  can  provide  the  same  graphic  displays  as 
well  as  distorted  shapes  of  the  structure  which  describe  the 
resulting  displacement  fields.  Detailed  information  about  these 
graphic  displays  are  given  in  section  3.6.  To  obtain  a  graphic 
display,  the  user  must  make  the  entry: 

RUN  PGRAPH 

Then,  PGRAPH  reads  the  JOBNAME.ATI  file  and  writes  down  three 
files  named  DJOBNAME.LST,  D JOBNAME . FOR  and  DJ0BNAME.COM.  Moreover, 
during  the  run,  questions  appear  on  the  screen  which  are  very 
similar  to  the  question  which  appears  during  a  PGEN  run.  The 
DJOBNAME.LST  file  is  identical  to  the  JOBNAME. LST  file  which  is 
obtained  after  a  PGEN  run  and  allows  a  detailed  check  of  the  data 
file  if  PGRAPH  is  run  at  the  preprocessing  level.  The  DJOBNAME.FOR 
file  contains  the  main  program  which  is  needed  to  obtain  the 
graphic  displays.  It  need  not  be  edited.  Finally,  the  DJ0BNAME.COM 
file  is  used  to  compile  and  link  the  program*.  Generally,  it  need 
not  be  edited. 

After  the  PGRAPH  run,  the  user  must  compile  the  main  program 
and  link  it  to  the  requested  libraries.  This  step  is  realised 
simply  using  the  entry: 

6D JOBNAME 

which  provides  the  executable  (DJOBNAME.EXE).  Finally,  the  user 
must  run  the  program  using  the  entry: 

RUN  DJOBNAME 

During  this  run,  the  code  asks  the  Job  name  (JOBNAME),  the  Job 
type  (ATI)  and  requires  the  choice  of  the  pre  or  postprocessing 
level.  The  preprocessing  level  can  be  used  as  soon  as  the 
JOBNAME.ATI  file  is  available.  The  postprocessing  level  can  only 
be  used  If  a  full  ATILA  run  has  previously  been  performed  (main 
run  or  restart)  and  if  the  necessary  storages  have  been  made.  The 
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whole  graphic  display  phase  must  be  performed  using  a  graphic 
terminal  and  the  Information  about  the  display  type,  its 
orientation,  its  title  ....  are  required  during  the  run  by 
questions  which  are  written  down  on  the  screen. 


*.  The  link  of  the  main  program  to  the  GKS  library  is  provided 
by  the  file  LNK1D.COM  which  is  used  by  DJ0BNAME.COM.  Before  the 
first  PGRAPH  run,  the  user  must  check  that  this  link  is  correct, 
taking  account  of  the  GKS  installation  on  the  corresponding 
computer . 
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5.5  THE  AUTOMATIC  MESH  GENERATION 


The  use  of  the  automatic  mesh  generator  MOSAIQUE  Is  described 
in  the  following  chapter.  It  requires  a  JOBNAME.MOS  file  which  is 
written  down  by  the  user.  Then  the  user  runs  MOSAIQUE  with  the 
entry: 

RUN  MOSAIQUE 

MOSAIQUE  provides  two  files  named  JOBNAME.LST  and  JOBNAME.ATI. 

JOBNAME.LST  contains  information  which  are  only  needed  If  a 
problem  appears  during  the  generation.  The  JOBNAME.ATI  file  is  the 
data  file  for  the  ATILA  main  run.  Depending  upon  the  contents  of 
the  JOBNAME.MOS  file,  it  can  be  used  immediately  or  can  be  edited 
by  the  user  to  add  the  last  entries  and  data,  generally  the 
loading  or  boundary  conditions. 
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Automatic  mesh 
generat Ion. 


Preprocessing. 


Main  run. 


5.6  SUMMARY 


r—  JOBNAME. MOS 

L»RUN  MOSAIQUE 

n 

JOBNAME.LST 

J. 

JOBNAME.ATI 

*j 

-i — JOBNAME.ATI 
U.RUN  PGRAPH 

D JOBNAME.LST 

«— 

r-   DJOBNAME.FOR 

4— 

U  DJ0BNAME.COM 

U.  SDJOBNAME  1 

i—  DJOBNAME.EXE*— ' 

-U>  RUN  DJOBNAME 

Interact Ive 

graphic  display 

-p  JOBNAME.ATI 
U  RUN  PGEN 

JOBNAME.LST 

,—  JOBNAME.FOR 

*_ 

U  J0BNAME.COM 

*-" 

U  eJOBNAME 

a 

p- JOBNAME . EXE 
-U  RUN  JOBNAME 

JOBNAME.LST 

JOBNAME. WRK 

JOBNAME. SY2 

JOBNAME. SY4 
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Restart 


JOBNAME.ATI 
JOBNAME.WRK 
J0BNAME.SY2 
J0BNAME.SY4 


*-*  RUN  PGEN 


E 


JOBNAME.LST  *— 
JOBNAME.FOR  *- 
J0BNAME.COM  *- 


©JOBNAME 

JOBNAME.EXE 
RUN  JOBNAME 


1 


JOBNAME.LST  #- 

JOBNAME.WRK  — 

JOBNAME. SY2  — 

JOBNAME. SY4  «- 


Postprocessing 


■JOBNAME.ATI 
JOBNAME. SY4 

RUN  PGRAPH  — 

D  JOBNAME.LST* 
D JOBNAME. FOR  *- 
DJ0BNAME.COM  *-» 
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®DJOBNAME 

DJOBNAME.EXE 

RUN  DJOBNAME 

Interactive  graphic  display 


Due  to  their  large  size,  the  files  .LST,  .WRK,  .SY2,  .SY4  are  not 
stored  as  successive  versions  If  they  are  provided  several  times 
with  the  same  job  name  but  are  automatically  deleted  by  the  code. 
If,  for  any  reason,  the  user  wishes  to  save  one  of  the  previous 
version  before  a  new  run  with  the  same  job  name,  it  has  to  rename 
the  corresponding  file. 
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CHAPTER  6 
AUTOMATIC  MESH  GENERATION 


THE  FINITE  ELEMENT  CODE  AT  I  LA  HAS  BEEN  DESIGNED  BY  THE  ACOUSTICS 
LABORATORY  OF  THE  "INSTITUT  SUPERIEUR  D* ELECTRON  I QUE  DU  NORD" 
(I SEN.  LILLE.  FRANCE)  FOR  THE  "GROUPE  D' ETUDES  ET  DE  RECHERCHE  EN 
DETECTION  SOUS-MARINE"  (GERDSM.  DCAN.  FRENCH  NAVY.  TOULON).  THE 
SOFTWARE  ENGINEERING  IS  PROVIDED  BY  THE  SINAPTEC  COMPANY  (LILLE." 
FRANCE). 
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6.1  INTRODUCTION 


As  mentioned  In  section  1.3.2  and  in  the  introduction  of  the 
third  chapter,  partial  regularity  of  the  mesh  can  often  allow  the 
use  of  an  automatic  generation,  which  provides  node  coordinates 
and  element  topologies.  In  the  ATILA  code,  this  work  is  performed 
by  a  specialized  preprocessor  named  MOSAIQUE.  Then,  the  user  has 
only  to  define  a  gross  splitting,  built  with  super  elements,  and 
to  select  their  automatic  splitting  into  finite  elements 
(Figure  1).  From  the  numerical  point  of  view,  each  super  element 
is  transformed  into  its  natural  axes,  using  the  isoparametric 
transformation.  Then,  it  is  split  into  finite  elements,  following 
the  corresponding  data  provided  by  the  user,  and  these  finite 
elements  are  transformed  back  to  the  global  axes  using  the  reverse 
isoparametric  transformation.  Finally,  redundant  nodes  are 
eliminated  by  searching  nodes  which  belong  to  a  same  given  small 
2D  or  3D  window. 

Entries  for  MOSAIQUE  are  provided  by  a  specific  input  file 
written  by  the  user  and  named  JOBNAME.MOS.  At  the  end  of  the 
MOSAIQUE  run,  the  output  file  which  contains  the  generated  node 
coordinates  and  topologies  is  named  JOBNAME.ATI  (see  chapter  5). 
Moreover,  entries  and  data  which  are  required  by  ATILA  but  are  not 
necessary  for  the  mesh  generation  can  be  also  read  by  MOSAIQUE,  if 
present  in  the  JOBNAME.MOS  file,  and  be  copied  in  the  JOBNAME.ATI 
output  file.  Thus,  starting  with  a  simple  input  file 
(JOBNAME.MOS),  MOSAIQUE  is  able  to  provide  the  whole  ATILA  data 
file  (JOBNAME.ATI).  Finally,  it  has  to  be  noted  that,  because  the 
number  of  nodes  which  have  to  be  defined  for  MOSAIQUE  is  often 
small,  this  procedure  allows  the  user  to  modify  very  easily 
geometrical  parameters  (length,  width  ...)  as  well  as  the  mesh 
steps. 

_Most  of  the  entries  and  data  contained  in  the  MOSAIQUE ,  data 
file  are  identical  to  the  ATILA  entries  and  data.  However,  the 
NODES  and  ELEMENTS  entries  are  modified,  while  two  new  entries  can 
be  used.  These  entries  and  data  are  described  in  section  6.2. 
Moreover,  the  splitting  of  a  super  element  into  finite  elements 
can  be  obtained  using  different  frames.  The  frames  which  are 
retained  in  MOSAIQUE  are  described  in  section  6.3. 
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FIGURE  1.  SCHEMATIC  DESCRIPTION  OF  THE  GENERATION  PROCESS. 
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6.2  MOSAIQUE  ENTRIES  AND  DATA 


6.2.1.  THE  HEADER. 

Information  contained  in  section  3.2  is  correct  for  the  header 
of  a  MOSAIQUE  data  file.  This  header  is  reproduced  as  the  header 
of  the  resulting  JOBNAME.ATI  file.  Thus,  it  will  be  the  header  for 
the  ATILA  job. 


6.2.2.  THE  ENTRIES. 

The  ATILA  entries  ANALYSIS,  CLASS,  NLOAD,  PRECISION,  LCPDDC, 
MATERIALS,  GEOMETRY,  GEOMETRY  POLARIZATION,  REDUCTIONS, 
FREQUENCIES,  RADIATION  and  PRINTING,  if  provided  in  the  MOSAIQUE 
data  file,  are  simply  read  and  copied  In  the  resulting  JOBNAME.ATI 
file.  Obviously,  they  are  not  used  by  the  mesh  generation 
procedure.  The  END  entry  has  to  be  provided  to  MOSAIQUE.  As  well 
as  for  ATILA,  it  defines  the  end  of  the  entry  list.  It  will  also 
be  copied  in  the  resulting  JOBNAME.ATI  file. 

The  NODES  entry  is  used  to  provide  the  coordinates  of  the  nodes 
which  define  the  super  elements.  Its  syntax  is  exactly  the  same  as 
for  the  ATILA  data  file  and  the  node  coordinates  have  to  be  listed 
in  the  node  number  ascending  order  associated  to  the  gross  mesh. 
After  the  MOSAIQUE  run,  this  NODES  entry  will  be  copied  In  the 
JOBNAME.ATI  file  with  the  coordinate  list  which  is  associated  with 
the  generated  nodes.  The  new  node  numbering  will  be  the  node 
numbering  provided  by  this  list.  It  has  to  be  noted  that  the 
numbering  of  the  nodes  belonging  to  the  initial  gross  mesh  is 
modified  during  the  generation.  Moreover,  as  well  as  for  ATILA, 
nodes  have  to  be  set  in  the  compulsory  order:  solid,  fluid, 
radiating.  Finally,  the  entry  NEWAXES  and  SCALE  can  be  used  to 
define  the  node  coordinates  for  the  gross  mesh  as  well  as  they  are 
used  for  ATILA. 

As  described  in  chapter  3,  the  ELEMENTS  entry  is  divided  in 
sets,  each  set  being  related  to  a  given  element  type*.  Successive 
sets  have  to  be  separated  by  a  dummy  line  and  the  whole  entry  must 
also  be  closed  by  a  dummy  line.  However,  with  respect  to  ATILA, 
the  content  of  a  set  is  modified  in  MOSAIQUE  by  replacing  each 
topology  line  of  the  type: 

N11  N12  N13  .  .  N1P 

by  a  set  of  lines  of  the  type: 

N11  N12  N13  .    .    .    .  N1P 

-1  ITYP  NK  NL  NM  NN  I 

1X1  IX2  IX3  .  .  .  .  IXNK 

IY1  IY2  IY3  .  .  .  IYNL 

IZ1  IZ2  IZ3 IZNM 
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where 


N11,   N12,   N13 N1P  are  the  actual   numbers  of   the 

corresponding  super  element  nodes, 

the  second  line  (-1  ITYP  NK  NL  NM  NN  I)  provides  compulsory 

information  about  the  splitting  of  this  super  element, 

ITYP  refers  to  a  given  frame  (see  section  6.3), 

NK  is  the  number  of  generated  elements  along  the  OX  natural  axis 

of  the  super  element, 

NL  is  the  number  of  generated  elements  along  the  OY  natural  axis 

of  the  super  element  (if  necessary), 

NM  is  the  number  of  generated  elements  along  the  OZ  natural  axis 

of  the  super  element  (if  necessary), 

NN  is  a  dummy  field,  kept  for  special  uses, 

I   is  a  control   integer  which  must  be  equal  to  0  if  the 

splittings  are  regular  and  1  if  at  least  one  splitting  is 

i  rregular , 

1X1,  IX2,  IX3,  ...,  IXNK  define  the  splitting  along  the  OX 

natural  axis  of  the  super  element. 

IY1,  IY2,  IY3,  ....  IYNL  and  IZ1,  IZ2,  IZ3 IZNM  provide 

the  same  Information,  respectively  for  the  OY  and  OZ  natural 
axes  of  the  super  element. 


Then,  along  the  OX  natural  axis  of  the  super  element,  the 
corresponding  side  is  divided  into  NK  different  parts,  the  lengths 
of  which  are  respectively  equal  to: 

IX1.L/IX,  IX2.L/IX,  IX3.L/IX,  ...,  IXNK. L/ IX 

where  L  is  the  total  length,  IX  is  given  by: 

IX  -  1X1  +  IX2  +  IX3  +  ...  +  IXNK 

and  1X1,  IX2,  IX3,  ...,IXNK  have  to  be  integers.  The  same 
splitting  rule  is  applied  for  the  OY  and  OZ  natural  directions.  If 
I  is  equal  to  0,  the  three  lines  which  are  devoted  to  the 
description  of  the  splitting  are  not  needed.  If  I  is  equal  to  1, 
these  three  lines  must  be  provided,  even  if  one  or  two  of  them 
correspond  to  a  regular  splitting. 

In  the  super  element  topology,  a  mid  side  node  can  always  be 
omitted  if  the  corresponding  side  is  straight.  Nevertheless,  if  a 
topology  line  Is  Incomplete,  it  has  to  be  completed  with  zeros,  at 
least  until  the  last  defined  node.  As  an  example,  If  the  mid- side 
nodes  12,  14  and  21  can  be  omitted  in  the  topology  line  1  3  5  9  12 
14  16  21,  this  line  can  be  written: 

12  5  9  0  0  16  0         or        1  2  5  9  0  0  16 

In  opposition  to  the  ATILA  case,  when  solid  and  fluid  elements  are 
used  together,  the  elements  must  be  provided  In  the  order:  solid, 
interface,  fluid,  damping.  This  order  can  be  modified  by  the  user 
in  the  JOBNAME.ATI  file  after  the  MOSAIQUE  run,  if  needed. 
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Finally,  two  new  entries  are  available  for  MOSAIQUE  with 

respect  to  AULA.  The  first  is  used  to  provide  information  related 

to  the  generation  of  the  damping  elements  and  must  be  written  in 
the  fol lowing  form: 

SORTFLU I D 
NFLU  NAC  NSH 

where  NFLU  is  the  actual  number  of  the  first  node  belonging  to  the 
boundary  surface  of  the  gross  mesh,  NAC  is  the  actual  number  of 
the  node  which  is  the  boundary  surface  centre  (acoustic  centre) 
and  NSH  is  a  dummy  field  kept  for  special  uses.  If  the  boundary 
surface  centre  does  not  belong  to  the  mesh  and  cannot  be 
identified  with  a  node,  NCA  has  to  be  equal  to  zero  and  the  origin 
of  the  global  axes  must  be  this  centre.  This  entry  is  compulsory 
if  damping  elements  are  used.  The  second  new  entry  is  used  to 
modify  the  size  of  the  2D  or  3D  window  which  allows  the 
elimination  of  redundant  nodes.  It  has  to  be  used  If  redundant 
nodes  are  still  present  in  the  JOBNAME.ATI  file  or  if  actual  nodes 
have  been  deleted.  It  has  to  be  written  under  the  form: 

MERGE 
XM 

where  XM  will  be  the  new  size  of  the  window  (XM  is  a  length  and 
must  be  provided  in  meters).  This  entry  has  only  to  be  used  if  the 
above  mentioned  problems  appear,  which  could  be  the  case  If  the 
mesh  is  very  Irregular.  Then  the  best  way  of  choslng  the  XM  value 
is  a  trial  and  error  method. 


In  fact,  the  element  type  is  related  simultaneously  to  the 
super  element  and  the  generated  finite  elements.  This  situation  is 
simple  when  the  super  element  and  the  generated  elements 
correspond  to  the  same  geometrical  shape  (splitting  of  a 
quadrilateral  super  element  Into  QUAD08E  finite  elements, 
splitting  of  a  triangular  super  element  into  AXIS06P  finite 
elements  ...)  but  can  lead  to  severe  difficulties  If  geometrical 
shapes  are  different  (splitting  of  a  quadrilateral  super  element 
into  TRIA06E  finite  elements,  as  an  example).  In  this  last  case, 
the  user  must  provide  an  element  type  which  corresponds  to  the 
geometrical  shape  of  the  super  element.  Then,  MOSAIQUE  writes  down 
this  wrong  element  type  In  the  JOBNAME.ATI  file  and  the  user  has 
to  modify  this  type  before  using  the  data  file.  Work  Is  In 
progress  to  avoid  this  problem  in  the  next  ATILA  version. 
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6.2.3.  THE  DATA. 

Data  can  be  simply  added  to  the  entries  In  the  JOBNAME.ATI 
file,  after  the  automatic  mesh  generation.  This  procedure  is 
classical,  since  the  node  numbers  which  are  necessary  to  write 
down  the  data  are  only  available  after  the  generation.  However,  in 
special  cases  (multiple  runs  with  the  same  topology  and  different 
dimensions,  for  example),  the  user  may  wish  to  write  down  these 
data  in  the  JOBNAME.MOS  file.  Then,  the  rules  are  the  same  as  for 
an  ATILA  data  file  and  the  corresponding  lines  are  simply  read  by 
MOSAIQUE  and  copied  into  the  JOBNAME.ATI  file. 
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6.3  THE  GENERATING  FRAMES 


The  splitting  of  a  super  element  In  simple  finite  elements  can 
be  realised  in  various  ways,  depending  upon  the  shape  of  this 
super  element  and  the  type  of  the  requested  finite  elements 
(quadrilateral  2D  elements  with  8  nodes,  20  node  hexaedras,  6  node 
ax isymmetr ical  interfaces  ...).  The  available  types  of  splitting, 
named  frames  in  the  previous  sections,  are  marked  by  the  ITYP 
control  integer  and  described  in  the  following  lines. 


ITYP  -  1 


A  linear  super  element  is  split  into  2  node  linear  finite  elements 
(SPRI02E)  following  the  diagram: 


.  0  ...  0  ...  0  . 

NK  elements  (NK  -  2) 


The  super  element  is  defined  by  nodes  1  and  2  and  the  control 
integers  NL  and  NM  are  dummy. 


ITYP 


A  linear  super  element  is  split  into  3  node  damping  elements 
(AXIS03M  or  AXIS03D)  following  the  diagram: 


.  0  ...  0  ...  0  . 
NK  elements  (NK  -  2) 


The  super  element  is  defined  at  least  by  the  nodes  1  and  2  but 
generally  also  by  its  mid  node,  due  to  the  curved  shape  of  the 
damping  boundary.  The  control  integers  NL  and  NM  are  dummy. 
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ITYP 


A  I inear 
(AXIS06I) 


super  element  Is  split 
following  the  diagram: 


Into  6  node  Interface  elements 


Fluid  domain 


...0...0...0.. 

.   0 

...0...0...0.. 

.    2 

So  1  i  d  dotna  1  n 

NK  elements  (NK  -  2) 

The  super  element  is  defined  by  the  nodes  1  and  2.  Its  mid  node  is 
required  only  If  the  interface  is  curved.  The  control  Integers  NL 
and  NM  are  dummy.  The  orientation  of  the  normal  to  the  interface 
from  the  fluid  to  the  solid  being  compulsory  for  the  generated 
interface  finite  elements,  the  orientation  of  the  super  element  is 
also  compulsory.  Nodes  1  and  2  are  defined  in  the  above  order  If 
the  Oz  natural  axis  is  the  outward  normal  to  this  paper  sheet. 
Then,  the  Ox  natural  axis  is  oriented  toward  the  left  (attention 
must  be  devoted  to  this  point,  since  this  numbering  rule  which 
applies  to  the  super  element  is  the  opposite  of  the  rule  which 
applies  to  the  finite  elements). 


ITYP 


A  quadrilateral  super  element  is  split  into  8  node  quadrilateral 
elements  (QUAD08E,  PLAT08E,  FACE08E,  AXIS08P,  QUAD08P,  AXIS08F  or 
QUAD08M)  following  the  diagram: 


NL 
elements 
(NL  -  2) 


0..0..0..0..I 
NK  elements  (NK  -  3) 


8 


Chapter  6.  Automatic  mesh  generation 


The  super  element  Is  defined  by  the  nodes  1,  2,  3  and  4.  Mid  side 
nodes  are  only  required  if  the  sides  are  curved.  The  control 
integer  NM  is  dummy. 


ITYP 


A  quadrilateral  super  element  is  split  Into  6  node  triangular 
elements  (TRIA06E,  PLAT06E,  AXIS06P,  TRIA06P,  AXIS06F  or  TRIA06M) 
following  the  diagram: 


3...   0   ...0...   0   ...4 


NL  .   . 

elements    0..0..0..0..0..0..0 
(NL  -  3)  .   . 


0...   0   ...0...0   ...0 


1..0..0..0..0..0 
NK  elements  (NK  -  3) 


The  super  element  is  defined  by  the  nodes  1,  2,  3  and  4.  Mid  side 
nodes  are  only  required  if  the  sides  are  curved.  The  control 
integer  NM  is  dummy.  This  case  is  related  to  the  comments  at  the 
end  of  section  6.2.2  and  care  must  be  paid  to  the  element  type  (as 
an  example,  the  element  type  can  be  QUAD08E  In  the  MOSAIQUE  data 
f i le  and  TRIA06E  in  the  ATILA  data  f  I  le). 


ITYP  -  6 


A  triangular  super  element  is  split  into  6  node  triangular 
elements  (TRIA06E,  PLAT06E,  AXIS06P,  TRIA06P,  AXIS06F  or  TRIA06M) 

following  the  diagram: 
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NL 
elements 
(NL  -  3) 


0     0 


0  .  .  0  .  .  0 


0     0     0     0 


0  .  .  0  .  .  0  .  .  0  .  .  0 


0     0 


0    0     0     0 


1..0..0..0..0..0..2 
NIC  elements  (NK  -  3) 


The  super  element  is  defined  by  the  nodes  1,  2  and  3.  Mid  side 
nodes  are  only  required  if  the  sides  are  curved.  The  control 
integer  NM  is  dummy. 


ITYP 


A  quadrilateral  super  element  is  split  into  16  node  quadrilatera 
interface  elements  (QUAD16I)  following  the  diagram: 


0    .    .   0    .    .   0   . 

.   0    . 
.   0    . 

.   0 
..4 

0.                       0. 
.0                       .0 

0. 
.0 

NL 
elements 
(NL  -  2) 

0. 
.0 

0.. 
.0   . 

•                                                        ■    • 

.   0   .    .   0. .    .   0    .    .0. 
.   0    .    .  .0    .    .   0   .    .  .0 
0. 
.0                       0. 
.0 

0.. 
1    . 

.  0   . 
.  0   . 

.   0.. 
.   0    . 

•                                                        *    • 

.   0    .    .0. 
.   0    .    .    2 

NK  elements  (NK  -  2) 
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The  super  element  is  defined  by  the  nodes  1,  2,  3  and  4.  Mid  side 
nodes  are  only  required  if  the  sides  are  curved.  The  control 
integer  NM  is  dummy.  Attention  must  be  paid  to  the  orientation  of 
this  super  element,  the  normal  to  the  interface  being  always 
oriented  from  the  fluid  toward  the  solid. 


ITYP  -  10 


An  hexaedrlc  super  element  Is  split  into  20  node  hexaedric  finite 
elements  (HEXA20E,  HEXA20P  or  HEXA20F)  following  the  diagram: 
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1..0..0..0..2 
NIC  elements  (NK  -  2) 


NM 
elements 
(NM  -  3) 


NL 
elements 
(NL  -  2) 


The  super  element  is  defined  by  the  nodes  1,  2,  3,  4,  5,  6,  7  and 
8.  Mid  side  nodes  are  only  required  if  the  sides  are  curved. 


1 1 


Chapter  6.  Automatic  mesh  ger 


ITYP  -  11 

A  prismatic  super  element  with  a  triangular  basis  is  split  Into  15 
nodes  prismatic  finite  elements  (PRIS15E,  PRIS15P  or  PRIS15F) 
following  the  diagram: 


NM 
elements 
(NM  -  3) 


NL 
elements 
(NL  -  2) 


1..0..0..0..2 
NIC  elements  (NK  -  2) 


The  super  element  is  defined  by  nodes  1,  2,  3,  4,  5  and  6.  Mid 
side  nodes  are  only  required  if  the  sides  are  curved. 


ITYP  -  12 


An  hexaedric  super  element  is  split  into  15  node  prismatic  finite 
elements  (PRIS15E,  PRIS15P  or  PRIS15F)  following  the  diagram: 
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Chapter  6.  Automatic  mesh  generation 
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NL 
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The  super  element  is  defined  by  the  nodes  1,  2,  3,  4,  5,  6,  7  and 
8.  Mid  side  nodes  are  only  required  if  the  sides  are  curved.  This 
case  is  related  to  the  comments  at  the  end  of  section  6.2.2  and 
care  must  be  paid  to  the  element  type  (as  an  example,  the  element 
type  can  be  HEXA20E  in  the  MOSAIQUE  data  file  and  PRIA15E  In  the 
ATILA  data  f I  le). 


ITYP 
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A  linear  super  element  is  split  into  3  node  ax isymmetr leal  shell 
finite  elements  (SHEL03E).  The  information  related  to  ITYP  -  2  is 
also  correct  for  this  case  but  the  numbering  of  the  generated 
nodes  and  the  topology  of  the  generated  elements  take  account  of 
the  specific  numbering  given  for  this  shell  in  section  4.3 
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